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Abstract. We present a model for spectral theory of families of selfadjoint operators, and their 
corresponding unitary one-parameter groups (acting in Hilbert space.) The models allow for a scale 
of complexity, indexed by the natural numbers N. For each n E N, we get families of selfadjoint 
operators indexed by: (i) the unitary matrix group U(n), and by (ii) a prescribed set of n non- 
overlapping intervals. Take O to be the complement in R of n fixed closed finite and disjoint 
intervals, and let L^{Vi) be the corresponding Hilbert space. Moreover, given B £ U{n), then both 
the lengths of the respective intervals, and the gaps between them, show up as spectral parameters 
in our corresponding spectral resolutions within L^{fl). Our models have two advantages: One, 
they encompass realistic features from quantum theory, from acoustic wave equations and their 
obstacle scattering; as well as from harmonic analysis. 

Secondly, each choice of the parameters in our models, n £ N, B £ U{n), and interval configura- 
tion, allows for explicit computations, and even for closed-form formulas: Computation of spectral 
resolutions, of generalized eigenfunctions in L^{iX) for the continuous part of spectrum, and for 
scattering coefficients. Our models further allow us to identify embedded point-spectrum (in the 
continuum), corresponding, for example, to bound-states in scattering, to trapped states, and to 
barriers in quantum scattering. The possibilities for the discrete atomic part of spectrum includes 
both periodic and non-periodic distributions. 
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1. Introduction 

The study of unitary one-parameter groups ([vN49]) is used in such areas as quantum mechanics 
([Bar49, Chull, AHMll] to mention a few), in PDE, and more generally in dynamical systems, and 
in harmonic analysis, see e.g., [DHJ09]. A unitary one-parameter group U{t) is a representation of 
the additive group of the real line M, t € M, with each unitary operator U{t) acting on a complex 
Hilbert space J^. By a theorem of Stone (see [Sto90, LP68, DS88] for details), we know that there 
is a bijective correspondence between: (i) strongly continuous unitary one-parameter groups U (t) 
acting on J(f; and (ii) selfadjoint operators P with dense domain in Jif. 

In quantum mechanics, the unit-norm vectors in the Hilbert space correspond to quantum 
states, and the unitary one-parameter groups U{t) will represent the solutions to a Schrodinger 
equation, with (Jff,U(t)) depending on the preparation of the quantum system at hand. In linear 
PDE theory, unitary one-parameter groups are used to represent time-dependent solutions when 
a conserved quantity can be found, for example for the acoustic wave equation, see [LP68]. In 
dynamical systems, selfadjoint operators and unitary one-parameter groups are the ingredients of 
Sturm-Liouville equations and boundary value problems. 

In these applications, the first question for {J^f, U{t)) relates to spectrum. We take the spectrum 
for U{t) to be the spectrum of its selfadjoint generator. Hence one is led to study {J^, U{t)) up to 
unitary equivalence. The gist of Lax-Phillips theory [LP68] is that {J^, U{t)), up to multiplicity, will 
be unitarily equivalent to the translation representation, i.e., to the group of translation operators 
acting in L^(IR, ^), the square-integrable functions from E into a complex Hilbert space The 
dimension of ^ is called multiplicity. For interesting questions one may take ^ to be of finite small 
dimension; see details below, and [JPTlla, JPTllb]. 

In this paper, we study a setting of scattering via translation representations in the sense of 
Lax-Phillips. 

To make concrete the geometric possibilities, we study here L^(f2) when is a fixed open subset 
of M with two unbounded connected components. For many questions, we may restrict to the case 
when there is only a finite number of bounded connected components in Q. 
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In other words, D, is the complement of a finite number of closed, bounded and disjoint intervals. 
We begin with Dirichlet boundary conditions for the derivative operator d/dx, i.e., defined on abso- 
lutely continuous functions with /' G and vanishing on the boundary of 17, / = on dil. 
Using deficiency index theory ([vN49, DS88]), we then arrive at all the skew-selfadjoint extensions, 
and the corresponding unitary one-parameter groups U{t) acting on L'^{^l). 

We expect that our present model will have relevance to other boundary value problems, for 
example in the study of second order operators, and regions in higher dimensions; see e.g., [Bra04]. 

1.1. Overview. In this setting, we resolve the possibilities for spectrum, and we show how they 
depend on the respective interval lengths, and their configuration, i.e., the length of the interval-gaps, 
as well as of the assigned boundary conditions. Our conclusions are computational, in closed-form 
representations; and expressed in terms of explicit direct integral formulas for each of the unitary 
one-parameter groups U{t). 

For each of these unitary one-parameter groups Usit), we compute its spectral decomposition as 
an explicit direct integral of generalized eigenfunctions ipx, functions of two variables, the spectral 
variable A S M, and of the spatial variable x G il. The direct dependence of generalized eigenfunc- 
tions on the boundary condition B is computed. The functions ^\ fall within the family known as 
exponential polynomials (see e.g., [AHDIO, MNIO]), or Fourier exponential polynomials. We further 
identify in detail those special selfadjoint extension operators for which there is embedded discrete 
spectrum. 

We now move on to the technical details in our construction, beginning with operator theory. 
A more detailed overview is postponed until the start of section 3. In fact, we will have a fuller 
discussion of applications in the three sections 3 through 5 below. In each, we begin with an outline 
of both the new main ideas introduced, as well as their spectral theoretic relevance to quantum 
mechanics, to wave equation scattering, and to harmonic analysis. 

In section 8, for comparison, we consider some cases when the give open set ft has an infinite 
number of connected components, still including the two infinite half-lines. This is of interest for 
a variety of reasons: One is recent studies of geometric analysis of Cantor sets [DJ07, DJll, JP98, 
PWOl]; so the infinite component case for J7 includes examples when fl is the complement in M of a 
Cantor set of a fixed fractal scaling dimension. This offers a framework for boundary value problems 
when the boundary is different from the more traditional choices. And finally, the case when the 
von Neumann-deficiency indices are ((xi,(xi) offers new challenges; see e.g., [DS88]; involving now 
reproducing kernels, and more refined spectral theory. Finally, these examples offer a contrast to 
the finite case; for example, for finitely many intervals (Theorem 3.20) we prove that the Beurling 
density of embedded point spectrum equals the total length of the finite intervals. By contrast, when 
ri has an infinite number of connected components, we show in sect 8 that there is the possibility of 
dense point spectrum. 

1.2. Unbounded Operators. We recall the following fundamental result of von Neumann on ex- 
tensions of Hermitian operators. 

In order to make precise our boundary conditions, we need a: 

Lemma 1.1. Let Q <Z M. be as above. Suppose f and f — -^f (distribution derivative) are both 
in L^(ri); then there is a continuous function f on fl (closure) such that f — f a.e. on J7, and 
limi^l^oo f{x) = 0. 
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Proof. Let p e M be a boundary point. Then for all a: G fJ, we have: 

f{x)~f{p)^rny)dy. (1.1) 
Jp 

Indeed, /' € L]^^ on account of the following Schwarz estimate 

\f{x)^ f{p)\<^/W^\\\f'\\^.^^y 

Since the RHS in (1.1) is well-defined, this serves to make the LHS also meaningful. Now set 

fix) :=/(p)+ r fiy)dy, 



and it can readily be checked that / satisfies the conclusions in the Lemma. ■ 

Lemma 1.2 (see e.g. [DS88]). Let L be a closed Hermitian operator with dense domain o, 
Hubert space. Set 

^± = {V'± e dom{L*) I L*V'± = ±#±} 
^(L) ^ {U : S)+ &^\U*U = Pc0^,UU* = Pcg_} (1.2) 

where Psi± denote the respective projections. Set 

^(L) ^{S\L<ZS,S*^ S}. 

Then there is a bijective correspondence between ^^(i) and S'iL), given as follows: 
If U Cz "^{L), and let Ljj be the restriction of L* to 

{^0 + /+ + Uf+ I ^0 e f+ € (1.3) 

Then Ljj G S'iL), and conversely every S G S'iL) has the form Ljj for some U G "^iL). With 
S G S'iL), take 

U:^iS-iI)iS + iI)-'\^^ (1.4) 

and note that 

(1) U G "^(L), and 

(2) 5 = 

Vectors f G domiL*) admit a unique decomposition f — y^o + f+ + /- where ipo G domiL), and 
f± G -for i/ie boundary-form B(-,-), we have 

zB(/,/) = (L*/, f)-{f,L*f) 

= II/+Il'-Il/-Il'- (1-5) 

1.3. Prior Literature. There are related investigations in the literature on spectrum and deficiency 
indices. For the case of indices (1,1), see for example [STIO, Marll]. For a study of odd-order 
operators, see [BH08]. Operators of even order in a single interval are studied in [Oro05]. The 
paper [BV05] studies matching interface conditions in connection with deficiency indices (m, m). 
Dirac operators are studied in [Sak97]. For the theory of selfadjoint extensions operators, and their 
spectra, see [Smu74, Gil72], for the theory; and [NazOS, VGT08, VasOT, Sad06, Mik04, Min04] for 
recent papers with applications. For applications to other problems in physics, see e.g., [AHMll, 
PR76, Bar49, MK08]. And [Chull] on the double-slit experiment. For related problems regarding 
spectral resolutions, but for fractal measures, see e.g., [DJ07, DHJ09, DJll]. 
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The study of deficiency indices (ri, n) has a number of additional ramifications in analysis: Included 
in this framework is Krein's analysis of Voltera operators and strings; and the determination of the 
spectrum of inhomogenous strings; see e.g., [DSOl, KN89, Kre70, Kre55]. 

Also included is their use in the study of de Branges spaces, see e.g., [Marll], where it is shown 
that any regular simple symmetric operator with deficiency indices (1, 1) is unitarily equivalent to the 
operator of multiplication in a reproducing kernel Hilbert space of functions on the real line with a 
sampling property Kramer). Further applications include signal processing, and de Branges-Rovnyak 
spaces: Characteristic functions of Hermitian symmetric operators apply to the cases unitarily equiv- 
alent to multiplication by the independent variable in a de Branges space of entire functions. 



2. Momentum Operators 

In this section we outline our model, and we list the parameters of the family of boundary value 
problems to be studied. We will need a technical lemma on reproducing kernels. 

By momentum operator P we mean the generator for the group of translations in L^(— oo,cxd), 
see (2.5) below. There are several reasons for taking a closer look at restrictions of the operator 
P. In our analysis, we study spectral theory determined by the complement of n bounded disjoint 
intervals, i.e., the union of n bounded component and two unbounded components (details below.) 
Our motivation derives from quantum theory, and from the study of spectral pairs in geometric 
analysis; see e.g., [DJ07], [Fug 74], [JP99], [LabOl], and [PWOl]. In our model, we examine how the 
spectral theory depends on both variations in the choice of the n intervals, as well as on variations 
in the von Neumann parameters. 

Granted that in many applications, one is faced with vastly more complicated data and operators; 
nonetheless, it is often the case that the more subtle situations will be unitarily equivalent to a 
suitable model involving P. This is reflected for example in the conclusion of the Stone-von Neumann 
uniqueness theorem: The Weyl relations for quantum systems with a finite number of degree of 
freedom are unitarily equivalent to the standard model with momentum and position operators P 
and Q. For details, see e.g., [J0r81]. 

2.1. The boundary form, spectrum, and the group U{n). Fix n > 2, let — oo < /3i < ai < 

/32 < Q!2 <•••</?„< q;„ < oo, and let 

Cn \ n 

\J[Pk,ak]\ Jk (2.1) 

fe=l / fc=0 

be the exterior domain, where 

Jo := (-oo,/3i), Ji := (ai, /32) , • ■ ■ , J„-i := (q!„_i, /?„), J„ := (a„,cx)). (2.2) 
Moreover, we set 

J_ := Jo, J+ := Jn (2-3) 
for the two unbounded components; see Figure 2.1 below. 
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J- - Jo J\ J2 

jSi cei P2 0-2 (63 



^3 Jn-l J+ - Jn 

/)4 a„.i j8„ a„ +00 



Figure 2.1. ft ^ UJJ^g Jj, = (^ULi -^fcj U (J_ U J+), i.e., 17 = the complement 
M of n finite and disjoint intervals. 



We shall write a — (oi) for all the left-hand side endpoints, and /3 — {f3i) for the right-hand side 
endpoints in dfl. 

Let L'^{n) be the Hilbert space with respect to the inner product 



n „ 

The maximal momentum operator is 

P 



1 d 



(2.4) 



(2.5) 



with domain !^{P) equal to the set of absolutely continuous functions on fl where both / and Pf 
are square-integrable. 

The boundary form associated with P is defined as the form 



z27rB(5,/) := {g,Pf)-{PgJ) 



(2.6) 



on S) (P). This is consistent with (1.5): If L = P^i„, then L* in (1.5) is P. Recall, ^(Pm»n) = {/ e 
^(P) ; / on dn}. 

Lemma 2.1. Lef cc = (ai), (3 = {l3i) be the system of interval endpoints in (2.2), and set 



Piif) :=/(/3) = 
for all f e 9{P); then 



V /(/?«) / 



, p2{f) :=/(a) = 



/ /(«i) \ 
/(a2) 



(2.7) 



z27rB(5,/) = (.g(a),/(a))c„ - (5(/3), /(/J))^. 
where (•, •)£„ is i/ie usual Hilbert-inner product in C". 

Proof. First note that for the domain of the operator L* in L'^{il), we have 

dom{L*) = {/ e . J/ ^ x2(f7)}. 

This means that every / € dom{L*) has a realization in C{fl), so continuous up to the boundary. As 
a result the following boundary analysis is justified by von Neumann's formula (1.5) in Lemma 1.2; 
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and valid for for all f,g e dom{L*): 

-i2nB{gJ) = {L*gJ),,-{g,L*f) 




n-l 

= (.g(/3),/(/3))c,.-(.g(a),/(a))c„. 

■ 

Corollary 2.2. It follows that the system (C",pi,p2), Pi{f) = /(/3) a^c? P2(/) = f{ct), represents 
a boundary triple, and we get all the selfadjoint extension operators for Pmin indexed by B € U{n); 
we shall write Pb- Explicitly, see e.g., [dO09], 

& [Pb] := {/ e {P) I Bp,{f) = p^if)} . (2.8) 

Remark 2.3 (The acoustic wave equation). Below we sketch the use of our interval-model for Lax- 
Phillips obstacle scattering ([LP68]) for the acoustic wave equation, with water waves; i.e., waves in 
a 2D medium. By [LP68], one knows that the solution to the wave equation, subject to obstacle 
scattering, may be presented by a unitary one parameter group U{t) acting on an energy Hilbert 
space J^E consisting of states representing initial waves as initial position and wave velocity. But, 
via a Radon transform (see [Hel98, LP68]), U{t), acting on the energy Hilbert space J^e, is in 
turn unitarily equivalent to a translation representation acting on L^(M,^). The Hilbert space ^ 
encodes the direction of the waves under consideration. In Figure 2.2 we illustrate a fixed compact 
planar obstacle, and four different states /i, /2, /a, and /4, each one with a different scattering profile. 
The first state /i transforms under U (t) in a manner unitarily equivalent to an interval model with 
two bounded component-intervals, see eqs (2.1) and (2.2). For the second state /2 the interval model 
has only one bounded component. The third state /a has no bounded component, but as with all four 
cases, the 17 model will have two unbounded infinite half-lines. The interval model for f^ corresponds 
to f2 = the complement in M of a single point. 

2.2. Some results in the paper. We identify a number of sub-classes within the family of all 
selfadjoint extensions Pb of the minimal operator in L^(r2). 

If the open set is chosen (as the complement of a fixed system consisting of n bounded, closed and 
disjoint intervals), then the set of all selfadjoint extensions is indexed by elements B in the matrix 
group U{n). The possibilities for the spectral resolution of a particular Pb are twofold: (i) pure 
Lebesgue spectrum with uniform multiplicity one; or (ii) still Lebesgue spectrum but with embedded 
point spectrum (within the continuum). 

While all the operators within class (i) are unitarily equivalent, it is still the case that, within each 
of the two sides in the rough subdivision, there is a rich variety of possibilities: Via a set of scattering 
poles, we show that the fine-structure of the spectral theory for each of the selfadjoint operators 
of Pb, and the corresponding unitary one-parameter groups UB{t), depends on all the geometric 
data: The number n, the choice of intervals, their respective lengths, and the location of the gaps; 
see Figure 2.1. More precisely, these spectral/scattering differences reflect themselves in detailed 
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/4 tangential case 



Figure 2.2. Obstacle scattering data as cross-sectional scans of a bounded planar object. 



properties of an associated system of scattering coefficients, see eq (3.1) in subsection 3.1 below.) To 
identifying particulars for a given unitary one-parameter group Usit) we study the location of a set 
of scattering poles. 

The resolution of these questions is closely related with a more coarse distinction: This has to 
do with decomposition properties for the unitary one-parameter groups UB{i) in L^(f2); a question 
taken up in the last three sections of the paper. 

In sections 2 and 3 below we prove the following theorem. 

Theorem 2.4. If B G U{n) is non-degenerate (see Definition 3.3), then there is a system of hounded 
generalized eigenfunctions {?/'^ ,A S M}, and a positive Borel function Fb{ ) on M such that the 
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unitary one-parameter group UB{t) in L'^{fl) generated by Pb has the form 

{UB{t)f){x)= f e^{~-t)U''\f) vf^(x)i^B(A)dA (2.9) 

for all f £ X £ Q, and t G M; where 

(^f\/) := / i/^)f{y)dy. 
^ ' Jn 

In section 2, we prepare with some technical lemmas; and in section 3 we compute explicit formulas 
for the expansion (2.9) above, and we discuss their physical significance. 

Our study of duality pairs x and A in systems of generalized eigenfunctions ipx is related to, 
but different from another part of spectral theory, dual variables for bispectral problems; see e.g., 
[Griill, GRIO, DG09]. 

Theorem 2.5. Let dasi-) measure in (2.9) and let Vb '■ iv^(O) L'^(R,<tb) be the spectral 

transform in (4-3) with adjoint operator Vg : i^(M, itb) ^ L'^{^). Then 

VbVb = and 
V'bVb = h^in)- 

Moreover, 

VBUB(t)V^ = Mt (2.10) 
where Alt is the unitary one-parameter group acting on L^(M, cts) as follows 

(Mtg) (A) = ex{-~t)g{X) 

for all <, A e M, and all g G i^(M, ctb). 

2.3. Reproducing Kernel Hilbert Space. In this section we introduce a certain reproducing ker- 
nel Hilbert space Jfi{il); a first order Sobolev space, hence the subscript 1. Its reproducing kernel is 
found (Lemma 2.6), and it serves two purposes: First, we show that each of the unbounded selfad- 
joint extension operators Pb, defined from (2.8) in sect 2.1, have their graphs naturally embedded in 
J^i{U,). Secondly, for each Pb, the reproducing kernel for J^i{U) helps us pin down the generalized 
eigenfunctions for Pb- The arguments for this are based in turn on Lemma 1.2 and the boundary 
form B from (2.7). 

Lemma 2.6. Let 

n 

n=\Jjk (2.11) 

be as above, and L'^(fl) be the Hilbert space of all L^ -functions on fl with inner product (•, and 
norm \\-\\q- Set 

J^iiQ) - {/ e L^n) \Df = fe L^{n)}:, 
then J^i(Q) is a reproducing kernel Hilbert space of functions on Q (closure). 

Proof. For the special case where = M, the details are in [J0r81]. For the case where 51 is the 
exterior domain from (2.11), we already noted (Lemma 1.1) that each / e J^i{ft) has a continuous 
representation /, and that / vanishes at ±oo. The inner product in Jlfi(fl) is 

{f,9).^,in)^{f:9)n + {f^9')n- (2-12) 
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Let X — UJJ^Q Jfe, and denote by J the interval containing x; and let p be a boundary point in J. 
Then an application of Cauchy-Schwarz yields 



fix) - f{p) = 25R / f(y)f'{y)dy 
Jp 

< 11/11' +ll/'ll'< 11/11^,(0) • 

We conclude that the linear functional 

J^i{n)3 f ^ fix) eC 

is continuous on J^ii^l) with respect to the norm from (2.12). By Riesz, applied to J^(f2), we 
conclude that there is a unique G J^(f2) such that 

/>) = (/c,,/),^^(^) (2.13) 

for all / e jfiin). 

If X in (2.13) is a boundary point, then the formula must be modified using instead fix^) — limit 
from the right if a; is a left-hand side end-point in J. If x is instead a right-hand side end-point in J, 
then use /(a;_) in formula (2.13). This concludes the proof of the Lemma. ■ 

Remark 2.7. If f2 is the union of a finite number of bounded components, and two unbounded, i.e., 
(See (2.1)-(2.3) and Figure 2.1. ) 

Vl = (-cx),^i) U y (aj,/3i+i) U (a„,oo) 



(2.14) 



where 
Set 



-OO < /3i < ai < /32 < "2 < • • • < "ri-l < Pn < ttn < OO. 



1 fc/3i 






( ka-^ 


\ 






k 




and k^ = 








} 









in ®"^]^ .yfiiVl). Let _B be a unitary complex n x n matrix, i.e., B G Uin); then there is a unique 
selfadjoint operator with dense domain l^iPs) in L^(il) such that 

r „ ] 

^iPs) = If eJfiin); f®- - j®l ±ikR-BkL) in 0J^^i(17) ^; (2.15) 
I. n times ' ^ 

and all the selfadjoint extensions of the minimal operator Dmin in L'^i^) arise this way. In particular, 
the deficiency indices are (n,n). 

Proposition 2.8 (Boundstates). Let n > 2; and set Ji = iai,f3i+i), J_ = (— oo,/3i), J+ = (q!„,oo) 
as in (2.14). Set fl = U^'^Ji, so 

L^iQ)^ L^in)(SL^i J-UJ+). (2.16) 
Of the selfadjoint extension operators Pb, indexed by B Cz Uin), we get the direct decomposition 

PB=Pn®Pe.t (2.17) 
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where P^^ is densely defined and s.a. in L'^(yt) and P^xt is densely defined and s.a. in L^(J- U J+), 
if and only if B (in U{n)) has the form 



( 



B 



\ 



(2.18) 



/ 



for some 9 G M/Z, and B G U{n - 1). 



Proof. Note that presentation (2.18) for some B e U{n) implies the boundary condition /(a„) — 
e{9)f{l3i) for / e !^{Pb) when Pb is the selfadjoint operator in i^(0) determined in Remark 2.7. 
And, moreover, the sum decomposition (2.17) will be satisfied. 

One checks that the converse holds as well. ■ 



Let B 



e U (n) , where u,w £ 



"^n — 1 



and c e C. In section 3, we consider the subset 



u B' 

c w* 

in U{n) given by u 7^ (see Corollary 3.28), but it is of interest to isolate the subfamily specified by 
u = 0. 

For 71 = 2, the unitary one-parameter group Usit), acting on L'^iVl), is unitarily equivalent to a 
direct sum of two one-parameter groups, Tp{t) and Tc{t). See Figure 2.3. These two one-parameter 
groups are obtained as follows: 

(i) Start with T{t), the usual one-parameter group of right-translation by / i— > /(• — t). The 
subscript p indicates periodic translation, i.e., translation by t modulo 1, and with a phase factor. 
Hence, Tp{t) accounts for the boundstates. 

(ii) By contrast, the one-parameter group Tc{t) is as follows: Glue the rightmost endpoint of 
the interval J_ starting at —00 to the leftmost endpoint in the interval J_|_ out to +00. These two 
finite end-points are merged onto a single point, say 0, on M (the whole real line.) This way, the 
one-parameter group Tc{t) becomes a summand of UB{t). Tc{t) is just translation in L^(M) modulo 
a phase factor e{ip) = e'^'^''' at a; = 0. 

For n > 2 (Fig 2.4), note the S-part {B e U{n — 1)) in the orthogonal splitting 



L\n) = L\\J'l-,'j,)(BL^m 

allows for a rich variety of inequivalent unitary one-parameter groups Ug{t). The case L^{Ji U J2) is 
covered in [JPTlla]. 

3. Spectral Theory 

In this section we establish a number of theorems giving detailed properties of each of the selfadjoint 
extension operators introduced in subsection 2.1 above. In Theorem 3.9 (the general case), we 
present the spectral resolutions as direct integrals: We give explicit formulas for the associated 
generalized eigenfunctions; and we study their properties. Among other things, we prove that they 
have meromorphic extensions to the complex plane C minus isolated poles, we give explicit formulas; 
and we study the scattering poles, both those falling on the real axis, as well as the complex poles. 

We now turn to some detailed spectral analysis of the operators acting on L^{^1). The first issue 
addressed may be summarized briefiy as follows: 
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barrier 1 


barrier 2 






f 7V(t) ' 







\ 1 




Ut), t > 



Figure 2.3. Infinite barriers (n = 2). Boundstates in one interval. 





barrier 1 




barrier 2 




barrier 3 




7_ 




Ji 




Ji 




^3 



/Si 



barrier n 



Jn-\ 




Figure 2.4. The complement of n bounded intervals in M (n > 2). Boundstates in 
the union of n — 1 intervals, and tunneling. 



We study three equivalent conditions 1 through 3 below, where: 

1. An clement B G U{n) is decomposable as a unitary matrix, i.e.. it has at least two non-trivial 
unitary summands Bi and B2. Note however, that this definition presupposes a choice of an ordered 
orthonormal basis (ONB) in C". 

2. As a self-adjoint operator in Pb is a corresponding orthogonal sum of the two operators 
Pi, i = l,2. 

3. The unitary one-parameter group Ub (t) generated by Pb decomposes as an orthogonal sum of 
two one-parameter groups with generators Pi, each unitary in a proper subspace in L^(f2). 

Some details about the corresponding summands in L^{^l), infinite vs finite. 
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The two infinite intervals: If a particular B in U{n) is decomposable, then the corresponding 
summands in (fi) arise from lumping together the spaces of the intervals Jj , j from to n, each 
corresponding to a closed subspace in L^(r2). But when lumping together these closed subspaces, 
there is the following restriction: one of the two infinite half-lines cannot occur alone: the two infinite 
half- lines must merge together. The reason is that for an infinite half-line, by itself yields deficiency 
indices (1,0) or (0,1). 

The finite intervals: If a subspace i^(Jj) for j from 1 to n — 1 occurs as a summand, there 
must be embedded point-spectrum (called boundstates in physics), embedded in the continuum. 

Caution about "matrix decomposition." The notion of decomposition for B in U{n) is basis- 
dependent in a strong sense: it depending on prescribing an ONB in C", as an ordered set, so 
depends on permutations of a chosen basis. Hence an analysis of an action of the permutation group 
Sn enters. So a particular property may hold before a permutation is applied, but not after. 

This means that some B in U{n) might be decomposable in some ordered ONB (in C") , but such 
a decomposition may not lead to an associated {Pb, L'^i^))- decomposition. 

For our matrix analysis we work with two separate notions, "non-degenerate" and "indecompos- 
able", but a direct comparison is not practical. The reason is that they naturally refer to different 
orderings of the canonical ONB in C". 

3.1. Spectrum and Eigenfunctions. Fix n > 2, let 51 be the exterior domain (2.1), see Figure 3.1 
below. 



Ao(A) 



J-=Jo 



A, (A) 



P2 



0-2 /?3 



A3(A) 



A„_,(A) 



«3 A a„-\ P„ 

Jl Jn-l 



AM 



J+=Jn 



FIGURE 3.1. ^ {^) = (ELo (A) XJ. i^)) ex (x) 



Let B = (bij) E U (n). Define the generalized eigenfunction by 



^ (^) E (^) ^''^ (^) (x) , A e 



\k=0 



where e\ (x) := e'^'^'^^^. The function 
given by 



a(-, •) : U{n) x M 
(S,A) := (4^' (A),...,4'') (A)) 



satisfies the boundary condition 

/ 4^^(A)eA(/3i) \ / 4^^(A)eA(ai) \ 



B 



4"^^ (A)eA (/32) 



V Ai^UA)eA(/3„) J 



(A)eA («2) 
V 4f ) (A) eA (a„) J 



with matrix-action on the LHS in (3.4). 



(3.1) 

(3.2) 
(3.3) 

(3.4) 
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Setting 

Dq (A) := diag (e\ (ai) , . . . , ca («„)) 
Dp (A) := diag (eA (^i) , . . . , ex (/?„)) 

and let 

(A) ■.= Dl{\)BDp (A) 
where i? is the matrix from (3.4). Then (3.4) can be written as 

/ (A) \ / ^1^) (A) \ 



Ba.fi (A) 



(A) 



V Af\ (A) / 



(A) 



V Ai^^ (A) y 

/ 4"' (A) \ 



(3.5) 



(3.6) 



where the matrix Ba^p{X) is acting on the column vector 



A[-^ (A) 



V (A) / 



In other words, with the 



definition (3.5), the two problems (3.4) and (3.6) are equivalent. 
Remark 3.1. Specifically, 

/ 6ii CA (/3i - ai) ••• 5i„eA (/3„ - ai) \ 
BcA>^)= ; ■•. ; (3.7) 

V bni e\ (/3i - a„) • • • 6„„ ex (^„ - «„) / 

3.2. The Role of U{n). The role of the group U{n) of all unitary complex matrices is as follows: 
On C" X C" (~ C^), we introduce the form S(-, •) from (2.7); 

B(z,C) = ||z|r-||Cf (3.8) 

where ||z||^ = is the usual Hilbert norm-squared. 

The projective space Pn,n is the complex manifold [Wel08] consisting of all complex subspaces 
L C C" X C" such that PriL = C", and 

B{z, C) = 0, for all (z, C) G L. (3.9) 

We use the notation Pri{zX) = z. 

The direction from U{n) to P„^„ is easy: If i? e U(n), set 

L{B):^{{z,Bz); zeC"}; (3.10) 

it is then clear that L{B) E Pn,n- 

For the converse argument, show that U{n) 3 B ^ L{B) maps onto P„,„, see for example [WelOS]. 



3.3. A Linear Algebra Problem. To understand the coefficients Ai{X) in the representation (3.3) 
of the generalized eigenfunctions, we will need a little complex geometry and linear algebra. 
Fix n > 2, and let 

u B' 

c w* 

where u,w E C"^^, and c E C. 



B 



E U{n) 



(3.11) 
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Definition 3.2. An element B E U (n) is said to be indecomposable iff it does not have a presentation 

Bi 



B = 



Bo 



(3.12) 



1 < fc < 71, Bi G U{k), B2 G Ui^n— fc); i.e., iff _B as a transformation in C" does not have a non-trivial 
spHtting Bi Q) B2 as a, sum of two unitaries. 

(The blank blocks in the block-matrix from (3.12) are understood to be a zero-operator between 
the respective subspaces. For more details, see section 6. ) 

Definition 3.3. Let B £ U{n) as in (3.11). We say B is degenerate if 1 e sp{B'), i.e., there exists 
C e C"-i\{0} such that B'C = C- 



Theorem 3.4. Let B 

the solution to 



u B' 

c w* 



G U{n) as in (3.11), where u,w E C" , and c G 



B 



Vl 



\ Vn-l J 



V2 



\vn J 



are as follows: 

(1) If B is non-degenerate: 



( ] 













{In-l-B'y^U 
W, {In-1 - B'y^ U 



for some constant Xg € C. 
(2) IfB IS degenerate: let C € ker(/„_i - B'), ( G C""i\{0}, then 
(a) If u not in the range of In-i — B' , 



(b) For u in the range o//„_i — B' , 




C 



( ] 













1 

Co 





c 



for some constant xq E C and some fixed (^q such that u (/n-i ~ -B')Co- 
Proof. Note (3.13) is equivalent to 



UVQ + B' 







( .1 \ 









and 



Then 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 
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/ ^1 



cvo + [w, 



(3.18) 



^ V Vn-1 

If 1 1^ sp{B'), solving (3.17) & (3.18) gives rise to (3.14). The remaining cases are similar. 
Example 3.5. Suppose n = 3, then 

10 
B = \ -1 

1 



(3.19) 



is degenerate. Here, B' — 



1 

-1 



Example 3.6. For n = 3, let 



and B'( — (, where C — 



0-1 
B = I -1 

1 



(3.20) 



so that £?' = ^ ^ 
Example 3.7. For n ~ 4, 



Note B'( = C, where C = 



1 



hence B is degenerate. 



B 



/ 1 \ 

1/2 1/2 

1/2 -1/V2 1/2 

V 1/V2 -1/V2 J 



is degenerate and 











h-B' 



1 



V2 2 
1 1 
\/2 2 



Hence u = | 1/2 | is in the range of I3 — B' and consequently we get an example for case (2)(b) 
1/2 

of Theorem 3.4. 
Example 3.8. For n = 2, let 



i.e., |a|^ + = 1. Suppose 



That is, 



B 



a b 
—b a 



a b 

—b a 



e SU{2), 



V2 



(3.21) 



avQ + bvi ~ Vl 
-bvQ + avi — V2 ■ 
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If 6 ^ 1 (non-degenerate), then 



vi \ = xo 

V2 



( ' \ 



a 
1-b 

1-b 



V H J 

If 6 = 1 (degenerate, a = 0) , the solution space is two dimensional, given by 



1 







Xo \ +2/0 1 , a;o,yo e C. 



-1 







3.4. The Generalized Eigenfunctions. We apply results in the previous section to the generalized 
eigenfunction in (3.1)-(5.1). 



Theorem 3.9. Fix B e U{n), and let 



(3.22) 



he the generalized eigenfunction in (3.1) satisfying the boundary condition (3.4)- Then a(i3,A) 
(^Aq^^ (A) , . . . , A^^"^ (A)^ in (3.3) is a solution to 



Ba,p (A) 



/ ^^"^ (A) \ 



V Af\ (A) / 



/ Af^ (A) \ 



V (A) ; 



(3.23) 



where B^^p — W^BD^, see (3.5) and (3.7). Moreover, writing 

«(A) B'^^piX) 



Ba.0 (A) = 



where 



c(A) w{X)* 
c(A) = 6„,ieA (/3i - a„) , 



«(A) 



and 



I hi &\ (/3i - ai) \ 
&21 Ca (/3i - 02) 

V ^1-1,1 ga (/3i - a„_i) / 
/ &I2 GA (/32 - ai) 



/ 6n,2 Ga (/32 - ««) \ 
&n,3 Ca (/33 - OLn) 

V ^ri.n Ga (^„ - a„) / 
&ln Ca (/3n - ai) 



"n— l,n GA (/3 n 



V ^ri-1,2 Ga (/3l - ftri-l) 

i/ien i/ie solution to (3.23) are as follows: 
Setting 

Ap = {A e M I det (/„_i - B'^^piX)) = O} 



(3.24) 



(3.25) 
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(1) If Ap — (p, then Ba^p{X) is non-degenerate, and 



I A^^' (A) \ 



V A!^n^ (A) ; 



{ln-1 - B'^^^iX)) U{X) 
^ c (A) + (w (A) , (/„_i - B'^^ (A)) « (A)^ ^ 



(3.26) 



for some constant Xq G C. The points A G Ap from (3.25) are the real poles in the functions 
Aj from (3.26). 

(2) Suppose Ap 7^ (f). For all X E Ap, Ba,p{X) is degenerate, and there is ( (X) e C"^^\{0}, such 



that C (A) e kcr (^/„_i - B'^ f^{X)j . Then 

(a) If u (A) is nof m the range o//„_i — ^(A) and (jn-i ~ B'^ Co (A) = u{X) 

I Af' (A) \ 



V A\f^ (A) / 





C(A) 
(t„(A),C(A)> 



(b) If u (A) is in the range of — B'^ ^{X) and 



I A^f' (A) \ 



V Ai"" (A) ; 



I Co (A) 

c{X) + {w{X)Xo (A)) 





C(A) 
(«;(A),C(A)) 



for some constant xq E C. In particular, Ap consists of eigenvalues for Pb. 

Proof. This follows directly from Theorem 3.4. ■ 

Corollary 3.10. Fix a system of interval endpoints a = (ai) and (3 = {Pi). Then the subset ofM. 

Ap = {A e M I dct (/„_i - B'^^piX)) = 0} (3.27) 
consists of isolated points, i.e., has no accumulation points. 
Proof. It follows from (3.24) that the function 

X^D{X)^ det - B'^^p (A)) (3.28) 

is entire analytic, i.e., is a restriction to M of an entire analytic function. 

To see this, note that A ^ B'^ ^(A) in (3.24) is entire; and since the determinant is multilinear, it 
follows D{-) in (3.28) is also entire. Since it is non-constant the properties of Ap (see (3.25)) follow 
from analytic function theory. ■ 

Corollary 3.11. Let fl be fixed as before, and select a B E U{n); then the functions Aj (•) in 
(3.26) and (5.1) have meromorphic extensions to C; the extension is obtained by replacing X in 
(3.24), (3.25) and (3.26) with z £ C The poles in the function C 3 z 1-^ A - (z) occur at the roots 



det (/„_i - B'^^piz)) = 



(3.29) 



and the embedded point- spectrum of the self adjoint operator Pg (in LF'{VL)) are the real solutions to 
(3.29). 
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Proof. The assertions in the corollary follow directly from the formulas (3.24) and (3.26) in Theorem 
3.9. ■ 

Remark 3.12. To find the meromorphic extension of the function 

M9 Ak^ -B;,^^(A))^' (3.30) 

from (3.26) in Theorem 3.9, we proceed as follow: Extend (3.30) by formally substituting z ^ C for 
A; and then proceed to compute the formal power series expansion for the function 

C3z^ i?„,^(z, B') (/„_i - D^i-z)B'Dp{z))-' (3.31) 

in the complement of the set of isolated poles. (The {n — 1) x [n — 1) matrix B' in (3.31) is fixed, 

but it is assumed to come from some i? = ^ ^ ^ E U (n) as in (3.11).) For iteration of the ^- 

derivatives in (3.31), it will be convenient to introduce = diag (oij)"^^, Lp — diag (/3j)J^2' 

and 

^cAM) MLp - L^M (3.32) 

defined for all (n — 1) x (n — 1) matrices AI. 

Then in the complement of the complex poles of Ra,p{z, B') in (3.31), we get 

-i d 
2-Ki dz 



R^^piz, B') - R^,p{z, B')5^^p{B')R^^p{z, B'). (3.33) 



And, as a result the higher order complex derivatives (j;;!^)" may be obtained from (3.33), and a 
recursion which we leave to the reader. It introduces a little combinatorics and an iteration of 5a, p 
in (3.32). 

In conclusion, we note that the complex extension 

C 9 z Ra,piz,B') 
is entire analytic in the complement of its isolated poles. 

Example 3.13. Let n = 2, and fix — oo < f3i < ai < /32 < a2 < oo. Let B = ^ ~ ) ' ^"'^^'^^ 
a,beC, and |a|^ + |6|^ = 1. Then 

D (A) = 1 - 6 exiP2 - ai), A e M. 

As a result, 

Ap = (f) |6| < 1 a^O. 
If a = 0, then there is a G M, such that b — e{9); and then 

Remark 3.14. Note that the complex poles discussed in Corollary 3.11 for Example 3.13 {b ^ 0) may 
be presented as follows: Select a branch of the complex logarithm "log"; then the complex poles are 

zeC,ze- ^——(^\ogb + z]\. (3.34) 

length (Ji) \2m ^ J j ^ ' 
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3.5. The Groups U{n) and U{n — 1). In the proof of Theorem 3.4, we considered the foUowmg 
operator/matrix block presentation of elements B G U{n), 



B = 



u B' 

c w* 



where u,w G C" ^, c G C, and B' is the (n — 1) x (n — 1) matrix in the NE corner in (3.35). 
We consider the coordinates in u as the matrix entries 

bii = Ui, 1 < i < n — \. 

For c G C, we have 

bni = c. 

The notation w* indicates that is a row-vector; we have 



(3.35) 

(3.36) 
(3.37) 



= Wj, 1 < j < n — 1. 



Finally we denote the Hilbert inner product (•, •) and it is taken to be linear in the second variable. 
With this convention we have 



Theorem 3.15. If B = 



u B' 

c w* 



ag{B) 



w u = {w, u) G C. 
G U {n) and g Cz U {n — 1), assume I ^ sp {B'). Then 
( gu gB'g-^ ^ 



c (gw) 



£U{n). 



(3.38) 



Uv = iv^)"^o G C"+i solve 



B 



then 















Vn-1 J 






( 




\ 



(3.39) 



V := 



V 



V2 



\ Vn-1 I 



(3.40) 



solves 



V2 



9 



V 



V Vn-l j 



( ( \ \ 

V2 



9 



\ Vn-1 j 
Vn 



(3.41) 
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Proof. Since B in (3.35) is in U{n), we get the following presentation of the C" norm: 

Wxqu + B'x\f + \cxo + {w, x)f = \xof + \\x\f 
^0 \ ^ j-.ri^ -y^g choose coordinates such that xq G C, and x £ C"^^. Since g e U{n — 1), 



for all 
we get 



for all 



Ixogu + gB'g ^x\\^ + {cxa + {gw,x)\'^ = \xo\^ + \\x\\'^ 



Xq 
x 



e C". The assertion in (3.38) follows from this. 



We now use Theorem 3.4 to solve the problem for ag{B). Hence the solution to the ag{B) 
problem is 





q 






/ 














= va 1 


V 


) 










) 





gu 

-1 



c + (gw, {In-i - gB'g 1) gu 



( 1 

g{In-i-B') 



which is the desired conclusion in (3.40). 

Inside the computation, we use the following formula from matrix theory 



and as a result 



- gB 
gw, (/„_i - gB'g-^ 



'"-y^gil^.^-Br'g-' 



gu 



gw,g{In-i - B') ^ u 



^(w,{In-i-B') u 



where we used g* g = In-i, i.e., g € U{n ~ 1). 



Corollary 3.16. Let B = 



u B' 

c w* 

diag{zjyy^ , Zj £ C, \zj\ < 1, then 



be such that, for some g € SU{n — 1), we have gB' g 



n-l 



det (/„_i - {gB'g-^)^^^ (A)) = \{{l z^ e{\Lk)) 



(3.42) 



k=l 



where Lj. = length{Jk), 1 < k < n. 

Lemma 3.17. Given B G U{n), then the following are equivalent: 
(1) ag{B) = B, for all g G U{n — 1), and 
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(2) B has the fo 



B 



In-l 

c 



, ce C,|c| = 1. 



Proof. Immediate from the definition of Ug, (3.38), i.e., 



an 



u B' 

c w* 



gu gB'g ^ 
c {gw)* 



(3.43) 



(3.44) 



see Corollary 3.10. ■ 
Lemma 3.18. // B is degenerate and ^ G C"^^ is an eigenvector of B' with eigenvalue 1, then 

PoC = PnC - 0. 

Proof Note B' = F^BP^ is contractive, and B'* ^ P^B'*P^^; and B'C, = ( implies that B'*C = C- 
Hence P,iC = C, P^C = C, and so F„C = ^bC = 0. ■ 

Corollary 3.19. Let B be degenerate. Then u and w are orthogonal to where is an eigenvector 
as above; i.e., 

{u,C)=0={wX)- 



Proof. By definition 

Bei = 

But recall P^C = PnC = 0, by Lemma 3.18. Then 

{u, = («, B'C) = {u, BC) = {B*u^ = (ei, = 
since PiC = 0. The same argmnent yields {w, C) = 0. 

Theorem 3.20. Set Jq — J^, J„ ~ J+. Let B be determined by (3.43) that is 

In-l 

c 



B 



, ce C, IcI = 1. 



(3.45) 



then the continuous spectrum of Pg is the real line and the discrete spectrum of Pg is [J^^l 
where £k = Pk+i — ctk is the length of the kth bounded interval. The multiplicity of each eigenvalue 
A is #{l<fc<n — 1| £fcA € Z} . Hence, is an eigenvalue with multiplicity n — 1 and counting 



multiplity the discrete spectrum has uniform density ^fe; sense that, for any a we hav 



number of eigenvalues in [a — n, a + n] 
2ri 



fc=i 



as n — > cxD. 
Proof. Note 



B 



/ Ao(A)eA(/3i) \ 
Ai (A)eA m 

V A„_i(A)eA(/3„) ) 



( A^{X)ex{ai) \ 
A2 (A) ex (02) 



\ A„ (A) ex (a„) / 



(3.46) 



is equivalent to 



Ak (A) ex {Pk+i) = Ak (A) ex (a^) , A: = 1, 
cAq (A) ex (/3i) = A„ (A) ex (a„) • 



, n — 1 
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Consequently, L'^{n) = L'^{J_ U J+) ® 0'^;^^ L'^{Jk) and 



where Pq is a selfadjoint operator acting in L^{J_ U J+) determined by c/(/?i) = f{an) and has 
Lebesgue spectrum, see Fig. 2.4, by [JPTlla], and P/j acting in L^i^Jk) is determined by f{Pk) = 
f{ak) and has spectrum -g^^^Z. Since the set has uniform density £, the density claim follows. ■ 



The same argument shows 



Corollary 3.21. // B 



u B' 

c w* 



where c = e{6i), B' = diag (e(02), • ■ • , then the con- 



tinuous spectrum of Pb is the real line and the discrete spectrum of Pb is ljfe=i 



^Z 



the multiplicity of each eigenvalue A is # {2 < fc < n | £kX — S/c G Z} , and counting multiplicities the 
discrete spectrum has density ^k- 

Remark 3.22. Recall the cyclic permutation matrix: 



S 



( 
1 



Vo 







1 0/ 



and S'^ =S* 



/O 1 

: 



Vl 



0\ 



1 

/ 



then 



BS 



u B' 
c w 



S 



B' u 

w* c 



(3.47) 



matrix product 

For an application of this Remark, see Section 5, cases 1 and 2 in subsections 5.1 and 5.2. 

u B' 

c w* 

conditions hold: 

\w\? P-w = In-1 
2 



Corollary 3.23. A n x n complex matrix 



is in U (n) if and only if the following list of 



B'*B' 

B'w + cu^Q 

||«;||' + |c|' = ll«ll' + ld' 



(3.48) 



1 



where the following notation is used for vectors x G C" ^. We denote the projection in C" ^ onto 
the one- dimensional subspace Cx by P^. 

Proof Combine (3.47) and (3.48). ■ 

Corollary 3.24. Let B G [/(n) have the representation given in Corollary 3.23 with entries, matrix 
corner B' , vectors u, w, and scalar c. Then B' is a normal matrix if and only if the vectors u and 
w are proportional, with the constant of proportion of modulus 1. 



Proof. Immediate from the system of equations (3.48). 
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Remark 3.25 (A dichotomy). Corollary 3.24 tells us precisely when B has its matrix corner B' a 
normal matrix. 

Combining Corollary 3.24 and Theorem 3.15, then note that, by the spectral theorem for normal 
matrices, we may pick g G U{n-1) in order to diagonalize the normal matrix B'; i.e., with gB'g^^ — 
diag(zi, ,Zn-i). 

For this case, we then get the following dichotomy: 

(i) The set Ap (= the real poles) is non-empty if and only the eigenvalue list {zj} contains an 
element of modulus 1. The corresponding selfadjoint operator Pb in L^(fi) then has embedded 
point-spectrum . 

(ii) If every number in the list {zj} has modulus strictly smaller than 1, then all the poles are off 
the real line, and as a result, Pb has purely continuous spectrum. 

Now the matrices in Corollary 3.24 account for only a sub-variety in U{n), but "large." Recall B', 
being a corner of a unitary, is typically not unitary; but, if it's normal, then the spectral theorem 
applies. 

3.6. Permutation Matrices. 
Example 3.26. Let 

/ 1 \ 
1 
10 

V 1 y 

That is, u = e^, w = ei, c — 0. Note B' is non-normal. We have 



and 



B = 



e U(4). (3.49) 




D{X) = dct [/3 - B'^^^iX)] = 1, VA e C, 



1 eA(/33-a2) ex{l33 + 134 - ai - a2) 
(h-B'^piX)) ' = ( 1 eA(/34-a2) 

1 



Hence, by Theorem 3.9 (1), we get 

Ai{X) = exiPi + (3^+ /34-ai-a2~ a^) 

A2{X) = ex{(3i+(34-a2-a3) (3.50) 
AaiX) = ex{(3i - aa) 

for all A G C. See Figure 3.2. 
Conclusions: 

(1) The functions A t-^ ^i(A), 1 < i < 3, have no poles in C. 

(2) The functions A i-> Ai{X), 1 < « < 3, are complex exponentials, depending only on interval 
endpoints. 

(3) The spectrum of Pb is purely continuous Lebesgue spectrum. 
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Figure 3.2. B is selfadjoint and B' is non-normal. 



Figure 3.2 illustrates the action of the unitary one-parameter group Usit) acting in L'^{Vl), and 
generated by the selfadjoint operator Pb coming from the boundary matrix B (€ C/(4)) from Example 
3.26, see eq (3.49)). 

The unitary one-parameter group Usit) in Example 3.26 acts by local translations to the right, 
acting on functions in fi. Action "locally" by translation here refers to translations to the right 
within the individual connected components in f2. Moreover, Figure 3.2 illustrates these local trans- 
lations when interval-endpoints are encountered. 

For comparison, we sketch, in Fig 3.3 below, the modification of the diagram (in Fig 3.2) when 
the boundary matrix B from eq (3.49) is changed into the 4x4 identity matrix I4. 

Conclusion: If i? = /4, then the associated group UB{t) is acting in _L^(il) by simply crossing over 
the gaps between successive components in f2, moving from the left to the right, jumping between 
neighboring boundary-points. 

Both the illustrations with the two versions of B correspond to a hit-and-run driver, constant 
speed, instantaneous jumps between components in f2. The second one {B = I4) rides right through 
without changing direction, but the first one is drunk and jumps between components in fl, in either 
direction, until eventually escaping to +00. 

Caution: With B — I4, the associated corner 3x3 matrix B' is still non-normal. In fact, as for 
the boundary matrix B in Example 3.26, the B' from B = I4 is nilpotent. 




Figure 3.3. B = I4, and B' is non-normal. 



Corollary 3.27. Let n > 2, fix as above; and let B G U{n) be a permutation matrix. Let Pb be 
the associated selfadjoint operator. Then the unitary one-parameter group UB{t) is acting in L^(i7) 
by local translations to the right of velocity 1. Let L be the sum of the lengths of the n — 1 bounded 
components in Q. Then, as t increases from —00 to +00, UB{t) acts in an interval of length L by 
simply crossing over the gaps between components in fl, jumping between boundary-points, in either 
direction. In time interval of length L, UB{t) makes a permutation of the n — 1 bounded components 
Jj in fl, where 1 < j < n. For every n > 2, there is a permutation matrix B G U{n) which makes 
the permutation of the intervals into the identity permutation, i.e., the local translations riding right 
through, jumps between successive components in D,, until eventually escaping to -\-oo. 
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Corollary 3.28. Suppose ^ " ^ G C/(n). Then 

B' (^U{n~l)^=^u = 0^w = 0. (3.51) 
Proof. Immediate from (3.48). ■ 

Corollary 3.29. Let n > 2, and let B e U{n). Consider the presentation B — (^^^ 
(3.35). Then the following bi-implication holds: 

c — Q B'* B' is a non-zero orthogonal projction in C"^^. (3.52) 
In particular, if (3.52) holds, then \\B'\\ = 1, and so B is degenerate. 
Proof. (=>) Assuming c = 0, then from the equations in the system (3.48), we get 

'l|w|| = lkll = i; 

B'*S' = /„_i-P^; (3.53) 

B'B'* = In-l — Pu, 

and so in particular, both B'*B' and B'B'* are orthogonal projections in C"^^. It is known that 
orthogonal projections have norm 1. so 

||B'||'-||i?'*B'|l = |l/„_i-P^|l=l 

as asserted. Indeed, the projection Pj; = — -Pm has rank n — 2 > 1 by the assumption in the 
corollary. 

The converse implication (<=) may be proved by the same reasoning. ■ 

Remark 3.30. We will show in section 5 that the unitary one-parameter group Usit) in i^(r2) has 
boundstates if and only if the condition (3.51) in Corollary 3.29 holds; see also Figure 2.4. 

Example 3.31. Let i7 = ^ ^ ) ^ SU{2), \a\^ + \b\^ = 1, z e C, \z\ = 1, and set B = 
orthogonal projection onto the one-dimensional subspace in spanned by ei. 

Lemma 3.32. Let n > 2, and let Ck — (a,) and (3 = {Pi) be a system of interval endpoints in 
^a,/3 = U"^QJi, where Jq = = (— oo,/3i), Ji — {ai,l3i), 1 < i < n, Jn — J+ — (a„,oo) and with 
interval length Li :— — ai, see Figure 3.4. Let B G U{n) have the form 



a 


6 





-b 


a 











z 



B 

then 



u B' 

c w* 



D{\):=det{I,,^i-B'^f,), (3.54) 
as a function on R (X CzR), only depends on the interval lengths Li, i — \,2, . . . ,n — 1. 
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■■ I3„ - a,, 



a„ 



Figure 3.4. Interval lengths. 



Proof. Write B' = (gij), for M 9 A ^ D{X) in (3.54), we then get the inside matrix as follows: 



gi2 ex{Li) 

922 e\{fi2 - OL2) 



913 exWa - "i) 
523 exiL2) 



\ 3n-U2eA(/32 - a„_i) 
As a result, for -D(A), we get: 

/ ^-9i2ex{Li) -giaexiPz- ai) 

-ff22 eA(^2 - 02) l-523eA(i2) 



gi,nexi(3n ~ ai) \ 

52,n-l exi(3n-l - 02) 
gn-l,nex{Ln-l) ) 

-gi.nSxiPn - ai) \ 
-ff2,n-l eA(/3n-l - ^2) 

1 - eA(-^-n-l) / 



(3.55) 



(3.56) 



V ~5n-l,2eA(/32-an-l) 

The conclusion now follows by induction: The determinant may be computed from its (n — 2) x (n — 2) 
sub-determinants, doing the computations entry- by-entry in the first row of the inside matrix in 
(3.56). ■ 

Remark 3.33. A second proof may be obtained from the following determinant identity: Let T be a 
k X k matrix, and let D.;, i = 1, 2, be two unitary k x k matrices; then the following formula holds: 

dct (/ - DITD2) = dct {DID2) dot {DiD*2 - T) . (3.57) 

The proof of the latter follows directly from a use of the multiplicative property of the determinant. 

From Figure 3.1, note that the numbers Gi — ai — /?; are the lengths of the gaps between the 
successive intervals, i.e., between Ji_i and Ji, I < i < n. Set Gtot — — total gap-length. For 

the unitary diagonal matrices Di and D2 in (3.57), take Di = -Da (A) and D2 = Dp{X), where 

/ e(Aai) ••• 

e(Aa2) ■■• : 



\ 







V 

/ e(A/32) 


V 





e(A/33) 



e(AQ;„_i) / 
\ 







e(A/3„) J 



and 



(3-a 



{\)^dmg([e{XL,)X;zl) 
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Then, for the determinant function M 9 A — > D{X) in (3.54), we get the following useful identity: 

D{X) = e{-XLtot) det [diag (e(ALj)) - B'] (3.58) 
where B' is the (n — 1) x (n — 1) matrix from (3.35). 

Corollary 3.34. The determinant D in (3.54) o, function only o/ A G M, (Li) G M"^"'^, i.e., 
D{\)=D{X,Li,...,Ln-i,B'). 

Corollary 3.35. Let 



3 D{X) ^ det diag{e{XLj))" 



B' 



(3.59) 



be the determinant factor on the RHS in (3.58). For j = 1,2, ... ,n—l, let Dj{X) he the determinant 
of the (n — 2) X (n — 2) sub-matrix obtained from 



diag {e{XLj))" ^ — B' 



by omission of its j row, and its j column. Then 

ji-i 

-DiX) = J2L,e{XLj)D,iX). 



1 d 
2TTi dX 



(3.60) 



(3.61) 



Proof. Differentiate (3.59), viewing the determinant as a multi-linear function on the 7i — 1 columns 
in (3.60). Applying j^j^ ^^'^ column in (3.60) yields the desired formula (3.61). To see this, 
note that the matrix in (3.60) is 



/ e(ALi)-6i2 -6i3 

-622 e(AL2) - 623 



-hin 



\ 



\ -bn-l. 



-bn-1. 



Corollary 3.36. Let ol = (a^) and (3 = be as above, i.e., the specified interval endpoints. Let 
B [n) be written as ( ^ ^* ) . Then the following two conditions are equivalent: 



(1) The (a, f3)-B problem is non-degenerate for all a and f3; and 



(2) IIB'II < 1 where \\-\\ is the 



^n—l 



operator norm. 



Proof. Recall that (1) is the assertion that M 9 A B'^ ^ (A) satisfies 

det (/„_i - B'^j, (A)) ^ (3.62) 
for all A e M. In other words, (1) states that for all A e M, 1 is not in the spectrum of B'^ ^ (A). 
Using (2) we see that this is equivalent to B'^ ^ (A) 



< 1. 



Corollary 3.37. Let B 



u B' 

c w* 



G U {n). Then w is an eigenvector for B'* B' with eigenvalue 



|c|^; and u is an eigenvector for B'B'* with eigenvalue |c|^. 
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Proof. From Corollary 3.23 we know that 

B'w + cu = B'*u + cw = 0. 

Now, apply B'* to the first, and B' to the second, the desired conclusion follows, i.e., we get the two 
eigenvalue equations: 



B'* B'w — |c| w; and 
B'B'*u^\c\^u. 



(3.63) 
(3.64) 



Corollary 3.38. Let B = 



u B' 

c w* 



G U (n); then 
WB'W > Id 



(3.65) 



Proof. The proof divides into two cases. First if it) = 0, then B' E U {n — 1) by Corollary 3.28, and 
so ll-B'll = 1, and (3.65) holds. Conversely, suppose w ^ 0; then by (3.63), |c| £ sp {B'* B'), but from 
operator theory, we know that 

IIB'II^ = niax{s e M; s e sp{B'*B')} ^ \\B'*B'f = \\B'B'*\f ; 

hence |c|^ < ||i?'||^ which is the desired conclusion (3.65). The inequality (3.65) may be sharp. ■ 

Lemma 3.39. Let n > 2, and let A be an [n — 1) x (n — 1) matrix, u,w E C"^""^, c E C Suppose 



B 



A 


u 


w 





E U{n). Then B ~ B* if and only if 



(1) A = A*, 

(2) u = w, 

(3) ceR, 

(4) ||nf + |c|^ = l, 

(5) A^ + \\uf P^ = In^i, and 

(6) u E ^{A + cIn-i). 

Proof. A direct computation. ■ 

Example 3.40. Consider the following selfadjoint unitary 3x3 matrix B and its cyclic permutation 
B where 

/ a b ' \ 



B = 



a + c 



where a, 6, c > 0, and g = Vl — c^. One verifies that if c is close to 0, then g is close to 1; hence b 
must be close to 0, and a closed to 1. But the norm of the corner matrix 



B' 
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is a {— its numerical range). Thus, the inequahty (3.65) may be strict. 

f u B' \ ~ / B' u \ 

Remark 3.41. Let B — [ ^ I and B = [ ^ I . If C is an eigenvector for B' with 



c w J \ w c 

eigenvalue e{0) for some real 0, then 

w* c y y w*c 

But B is unitary, in particular ^ q ^ ^^^^ ^ ^w*^ ) have the same norm, hence 

w*C = 0. (3.66) 

Slightly generalizing a claim in Corollary 3.19. If B is selfadjoint, then w = u, hence implies (3.48) 
B'u = — cu. In particular, u is in the range of — B' , if c 7^ —1. On the other hand, if c = — 1, 
then (3.66) implies u = 0. 

Consequently, if B is selfadjoint, then we are never in case (2) (a) of Theorem 3.4. 

4. The Continuous Spectrum Is Simple 

The generalized eigenfunctions studied in the previous section (see Theorem 2.4 and eq. (3.1)) 
yield separation of variables, a harmonic part (in the spatial variable x as e\{x)), and a finite family 
of scattering coefficients {Aj (A)}, functions of the spectral variable A. In this section we study the 
meromorphic extension of scattering coefficients; extension to non-real values of A. 

We show (Theorem 4.1) that, if the first of the scattering coefficients is normalized to 1, then the 
continuous part of the spectrum for each of the operators is purely Lebesgue, with spectral measure 
having Radon-Nikodym derivative equal to the constant 1. We further show that each point on the 
real line M occurs in the continuous spectrum with multiplicity 1. 

Let the open set 57 be as before, i.e., n is the complement of n bounded closed and disjoint 
intervals. The minimal momentum operator will then have deficiency indices (n, n), and as a result, 
the boundary conditions are indexed by the matrix group U{n). As before, we denote the unbounded 
selfadjoint extension operators Pb indexed by a fixed element B € U{n), and the corresponding 
unitary one-parameter group is UB{t). These operators are acting in the Hilbert space L'^{n). 

We restrict the element B as in Theorem 3.9, i.e., it is assumed non-degenerate. In this generality 
we are able to establish (Thm 4.1) the complete and detailed spectral resolution for Pg, and therefore 
for the one-parameter group Usit) as it acts on the Hilbert space L?{^). We show that, if the first 
coefficient in the formula for the generalized eigenfunction system in (5.1), is chosen to be 1, then 
the measure ii^ the spectral resolution for UB{t) becomes Lebesgue measure. Moreover, we show 
that the multiplicity is uniformly one. In the theorem, we further compute all the details, closed 
formulas, for the spectral theory. 

Theorem 4.1. Let a = (a,;) and (3 = (/?.;) be a system of interval endpoints: 

" (X < f3i < ai < (32 < ■ ■ ■ < l3n < an < 00, (4.1) 

with Jo = J- = (~co, /3i), J„ — J-f- — {an, 00), and Ji = (a;, A+i), i = 1, . . . , n — 1. Let B E U (n) 
be chosen non-degenerate (fixed), and let 

(x) ) {x) = (j2 {x) (A) ) eA {x) (4.2) 
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be as in The 



are chosen as in (3.26) with Aq 
For f £ L"^ {ft), setting 



3.9, where n = \J"^^ J„ 

Xi XJi> ^ i ^ n, and where the functions f A| (•) j 



{VBf)iX) = {i'xJ)n^ J i^x{y)f{y)dy, 
we then get the following orthogonal expansions: 

/= / {VBf)Wi^x{-)d\ 



where the convergence in (4.4) to he taken in the Li^ -sense via 



Moreover, we have 



whe 



(o) = / \{VBf)(\)Vd\ feL'i^). 

JR 

VBUBit) = MtVB, teR 
(Mtg) {X)^ex{-t)g{X) 



for all t,XeR, and all g e L'^ (M) . 



(4.3) 
(4.4) 

(4.5) 
(4.6) 



i2 ^n)^L^ (n) 



Vb 



Vb 



(M) ^ L2 (M) 

Figure 4.1. Intertwining 

Remark 4.2. The reason for the word "generaUzed" referring to the family (4.2) of generahzed 
eigenfunctions is that, for a fixed value of the spectral parameter A, the function is not in 
L^(f2), so strictly speaking it is not an eigenfunction for the unbounded selfadjoint operator Pb 
in L^(r2). But there is a fairly standard way around the difficulty, involving distributions, see e.g., 
[JPTllb, Mau68, Mik04]. 



Example 4.3. Set n = 2, B 



_ j, a, 6 G C, |a| + |6| = 1. With the normalization 

j4q = 1, we get the following representation of the two function M 9 A i— (X), i = 1,2: Fix 
— oo < /?! < tti < < a2 < oo; set L := P2 — cti, and G := a2 — Pi, then 



A'^\X)^ 
4"Ha) = 



, and 



l-bex{L) 

ex{L-G)-bexiG) 
l-bex{L) 



(4.7) 
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Note the poles in the presentation of the two functions in (4.8). In the meromorphic extensions of 
the two functions, we have, for z e C, 

(S), . _ ae(z(/3i - aQ) 
^1 1-beizL) 

4.8 

AB), . _ e{z{L~G))-he{-zG) 
^2 ^'^^ l-he(zL) • 

Remark 4.4. It follows from Corollaries 3.10 and 3.11 that also in the general case with Vl open 
and associated ol = {cti)'^^i, and (3 = (/3i)"=ii the scattering coefficients A)^ (•) have meromorphic 
extensions. With this information, one may derive a suitable de Branges Hilbert space of meromorphic 
functions [dB84, ADV09], for a detailed and geometric analysis of the general case. It follows that 
formulas (4.7)-(4.8) in Example 4.3 are indicative for the study of the general case; only in the 
general case n > 2, the extension of the meromorphic functions C 9 z i-> (z) is substantially 
more difficult. 

Proof. Outline of proof in sketch. Given B in U{n), we get a specific selfadjoint operator Pb, as 
outlined in section 1 . And there is therefore an associated strongly continuous unitary one-parameter 
group Usit) generated by Pb and acting on i^(il). We begin with an application of the abstract 
spectral theorem: Given the selfadjointness of the operator Pg, we may apply the spectral theorem 
to it, but this yields only the abstract form of the spectral resolution; not revealing very much 
specific information: At the outset, the general theory does not say what the spectral data are, 
such as detailed information about the measure ub arising in the direct integral representation for 
Pb- Given the properties of Pb it does say that cjb must be absolutely continuous with respect 
to Lebesgue measure on M. But it does not say what the Radon-Nikodym derivative Fb is. Our 
assertion is that with the normalization Aq — 1, ■we obtain Fb — 1. 

Having ctb, we proceed to apply the theory of Lax-Phillips [LP68] to the unitary one-parameter 
group UB{t) as it is acting on L^(J7). To do this, we must assume that B is non-degenerate, so 
that Pb will have no point-spectrum. To apply Lax-Phillips, we do not need to know details about 
the measure cjb- Its abstract properties are enough. To begin with, we first establish that L'^{Jq) 
serves as an incoming subspace Si_ in L?{^) for the action of unitary one-parameter group UB{t). 
Recall Jo — (— oo, /3i). To show that this incoming subspace does satisfy the Lax-Phillips axioms, 
again we do not need detailed information about the measure ctb. Finally, with an application of 
Lax-Phillips, and a number of other steps, in the end, assuming Aq = 1, we are able to conclude that 
the measure ctb is Lebesgue measure on M. 

Proof in detail. First note that the assumption on B rules out point spectrum; see Corollary 
3.10. Using [JPTllb] and the theory of generahzed eigenfunctions [DS88, Mau68, Mik04, MM63], 
we note that there is a Borel measure ctb (dX) on M, absolutely continuous with respect to Lebesgue 
measure such that the formulas (4.4) and (4.5) hold with aB {dX) on the RHS of the equations. Our 
assertion is that aB {dX) = dX = Lebesgue measure on M, i.e., then the Radon-Nikodym derivative 

daB (dX) 



dX 



Fb (A) = 1. (4.9) 



The validity of (4.9) uses the assumption Aq = 1 in an essential way. Hence in (4.2), we have 

V'f ) (x) = I X(-oo,/^,) {x) + J2 (A) X, {x) 1 (x) . (4.10) 
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We will now suppress the i?-dependence in tp)^ (•) and (•). It is understood that ip\{-) is a 
function on fl, and each Aj (•) is a function on M; see Theorem 3.9 for the exphcit formulas. 
In the computation below, we will be using the normalized Fourier transform and its inverse 
Let Pj = multiplication by Xj for ^ j ^ '^j viewed as projection operators in (fi). We then 
have 

n 

Y,Pj=I = h^n). and P.P,. = <5,- ^P^-. (4.11) 

3=0 

From (4.10), we then get the following expression for Vb '■ (il) — > (cs): 

n 

{VbI) (A) = (Po/)^ (A) + J2M^){Pjff (A) , (4.12) 

i=i 

for all / e L2 {il), and all A G M; and 

/ = ^o/ + E^.(-)(^(^./r) ■ (4-13) 



(It is understood in (4.12), (4.13) and the sequel that Aj ~ Aj ' depends on a choice of B e U{n).) 

With B E U (n) specified as in the theorem, we get a unique selfadjoint operator Pg in (ft) 
as a selfadjoint extension of the minimal operator i.e., the minimal operator specified by the 

condition: / G (f]), /' e (Q) and / = on afi,'see [JPTllb]. 

Let, for t e M, 

Ub {t) : (f^) (1^) (4.14) 

be the corresponding strongly continuous unitary one-parameter group generated by Pg; see [Sto90, 
vN32, vN49]. 

Applying (4.13), we get the following formula for Ub (t) /, / G (il), understood in the sense of 
L^-convergence: 

(Ub (t) /) (x) = xo (x) (Po/) (x - i) + E (^) (^^- (•) (P^f)"")^ - ^) (4-15) 

for all / G and all cc G fi, i G M. 

We now prove (4.6) as an operator-identity, i.e., the assertion that that 

Vb : (n) ^ (as) (4.16) 

intertwines the two unitary one-parameter groups specified in (4.6). 
Let f £ L'^ (n), X e n, and A, < G M. Then, 

{VbUb (t) /) (A) 

= (V^A,(7s(i)/)o (by (4.3)) 

= {Ub {—t) V-'A, /)o (by the theory of generalized eigenfunctions) 

= (eA(t)^A,/)^, (by (4.10)) 

= {~t) {iPxt f)fi (since (•, •) is conjugate linear in first variable) 

= ex{-t){VBf){X) (by (4.3)) 
= iMtVBf){\). 
Since this holds for all A G M, the desired formula (4.6) is verified. 
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We now establish formulas (4.4) and (4.5) first for f G (Jo) = (— oo,/?i); and we recall from 
[JPTllb] that this subspace serves as an incoming subspace for U b (t) in the sense of Lax-Phillips 
[LP68]; see also [JPTllb], i.e., 



^_ - ^^0 = (Jo) = (-oo,/3i) . 
Proof of (4.4): For /o € and x E il; then (in the sense of i^-convergence) : 
/o (x) = XJo (x) fo (x) 

= XJo (x) f (x) (Po/o)^ (A) d\ (by (4.12)) 



(4.17) 



XJo {x) / {x) (Po/o)" (A) dX (by (4.10)) 

JWi 

{VbIo) (A) Va (x) dX (by (4.3) and (4.10)) 
which is the desired formula (4.4). 

Proof of (4.5): By the spectral theorem (see [JPTllb]), the measure (Jb (A) satisfies 

Wflllnn) = f \{VBf)iX)\'c7B{dX), 



(4.18) 



see also (4.16). Now specialize to f ~ Jq E C (fi). Using (4.4), and Parseval's formula, we get 



l/ol 



L2(r!) 



dX 



(4.19) 



fo (A) aB (dX) / fo (A) 
which is (4.5) on vectors fo G To see this, use (4.12). 

= may be stated in terms of the projection- 

Eb (•) : {Borel-sets in M} {Projections in (Q)} 



The conclusions (4.4) and (4.5) for vector fo e 
valued measure 



as follows: 



\\EB{dX)fo\\l^ /o(A) 



dX. 



It remains to prove that 

\\Eb idX)f\\l = \{VBf)iX)fdX 

holds, for all / e L^{n). 

But by Lax-Phillips [LP68[ and [JPTllb[, the linear span of the vectors 

{UB{t)fo;teRJoe&-} 



(4.20) 
(4.21) 

(4.22) 



is dense in {fl). This is where non-degeneracy of B is used. 

As a result, it is easy to establish (4.21) when f E (fi) has the form f = Ub (t) fo, t £ 

foe^- = J%. 

We proceed to do this. We have: 



\EBidX)UBit)fo\\l = \\EB{dX)fo\\l: 



(4.23) 
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and for the RHS in (4.21) with f = Ub {t) fo- 

\{VbUb (t) fo) {X)\^ dX 
= \exH){VBfo){X)\'dX (by (4.6)) 

fo{X)fdX; (by (4.12)) (4.24) 

The two right-hand-sides in the last two equations (4.23) and (4.24) agree as a consequence of (4.20), 
and we can therefore conclude that (4.21) holds for all / = Ub {t) fo as asserted. ■ 

Remark 4.5. The axioms for Eb (•) in (4.20) and (4.21) are as follows: 

(i) Eb (S) is a projection in (17) for all Borel subsets 5 C K, 5* G ^; 

(ii) S§ 3 S ^ Eb (S) is countably additive; 

(iii) Eb (S, n ^2) - Eb (S,) Eb (82), V^i, ^2 € ^; 

(iv) f = J^ Eb (dX) f holds for all f £ {n): 

(v) Ub (t) / = Jj^ex {-t) Eb (dX) f holds for all f e (fl), t <eR. 
The conclusion in (4.21) may be restated as follows: 

For Borel sets 5 (e ^), let Ms := multipHcation by xs in L'^ i^), and let Vb : (fl) ^-^-^ 
be the transform in (4.3) and (4.4); then 

EbIS) = V^MsVb; Se£§. (4.25) 

Convention. For functions / on E, we set Ma to be the corresponding multiplication operator 
(Maq) (A) = A {X) g (A), A S K, with adjoint M\ = , and ~ denoting complex conjugation. On 
L'^{ft) C i^(K), we view the Fourier transform as a unitary operator so J^f — f, and = g^, for 

/ B) 

Corollary 4.6. Let fl, and B G U{n) be specified as in the theorem, and let {V'^ } be the system 

( B') 

of GEFs in (4-10) with coefficients {Al }"=o. Then the spectral transforms Vb and Vg from (4-3) 
have the following representation as Fourier integral operators: 



Vb = ^,9^*M\^.^Pj, and (4.26) 

n 

V2 = ^P,=^*M^^.^ (4.27) 



where Pj = M^j , < j < n, and recall 
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Proof. For all / G L'^i^), by (4.3), we have 



(^b/)(A)= / ^Px{y)f{y)dy 



J2x^ {y)Af^ (A) {y)fiy)dy 

\i=0 ) 

^YA^)l ^)x^ {y)f{y)dy 



= Y.Af'> (A)^(P./); 
and this yields (4.26). On the other hand, for all g e L^(M), we have 



l\J2x^ ix)Af^ (A)j e, ix)gi\)d\ 
^X.(^) / AP^(A)e,(x)5(A)dA 



i=0 

which gives (4.27). 

Remark 4.7. For relevant details on Fourier integral operators, see e.g., [Duill] 
Corollary 4.8. Select a pair of elements B and C in U{n) specified as in Corollary 4.6. Let {A\ 



(B) 



and {A^- ') he the corresponding systems of scattering coefficients. Then for the operator VqVb 



have 



VSVb = P^■^*APAf^^Py, (4.28) 

i.e. in each (i^ j)- scattering block VqVb has the function 



Af\\)Af\\) (4.29) 

as a Fourier-multiplier. 

Proof. Immediate from (4.26) and (4.27) in Corollary 4.6. ■ 

Corollary 4.9. Fix 1 <i < n—\, and let Pi he the projection from L?{H) onto L?{Ji) — (a^, ft+i). 
Then for all f G L? (fJ), we have 

(A) - / Shann, (A - |A,|' (0 {P^ff (0 d^, (4.30) 

where 



Shanni{S^) := j {x)dx 

a.1 + Bi+i \ sin (7r£ (ft + 1 - a,)) 
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is the Shannon kernel on the bounded interval Ji; see [DM72]. 
Proof. By (4.12), we have 




(P,/) {x) =PJ (VbPJ) (A) (x) dX 



= X^ i^) f (VbPJ) (A) {x^ (x) i^x (x)) dX 



= X^{x) [ \A{X)\^{PJ)''{X)ex{x)dX. 




(P,/)^(A)= / X^{X~0\A^m'iP^ffiOd^ 
JR 



and (4.30) holds. 



4.1. An Inner Product on the System of Boundary Conditions. While large families within 
the selfadjoint extensions Pb, B £ U{n) are unitarily equivalent, there are much more refined mea- 
sures that pick out specific scattering theoretic properties for the selfadjoint operators and the cor- 
responding family of unitary one-parameter groups. Below we compute two such; one is an inner 
product, or a correlation function, defined initially on U{n) and then extended by sesqui-linearity. 
The second is the family of scattering semigroups, see [LP68]. 

Corollary 4.10 (An Inner Product on the System of Boundary Conditions.). Let a = (ai)f^i and 
(3 — {l3i)f^i be fixed as above. For each of the finite intervals Jj :— (a^, /3j+i), j = 1, . . . ,n — 1 in fl, 
let Shj — Shj. be the corresponding Shannon kernel 



^ ' JR 

For elements B G U{n) specified as in Theorem ^.1, let Vb and Vg be the associated transforms. 
For pairs of elements B,C Cz U{n) we have 




(4.32) 



For a pair of bounded Borel functions gi and g2 on M, set 




(4.33) 




(4.34) 



and that 



n-l 




(4.35) 



where Xfin '■= Y^^Zi Xjj = indicator function of the union U^Zi^j of the finite intervals. 
Note that (4-35) extends by sesqui-linearity to a Hilbert inner product {B,C); and then 




3 = 1 



Proof. Follow from Theorem 4.1 and Corollaries 4.9 and 4.9. 
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Example 4.11. Let n — 2. Fix a system of interval endpoints 

—oo < /3i < ai < < ^2 < oo 
and let = Jq = (— oo,/3i), Ji = (ai,/?2), and J2 = Joo = (0^2,00). 



Let B = 



a b 
—b a 
By (4.7), we have 



and C 



c d 

—d c 



where a,b,c,d G C, and |a|^ + = |c|^ + \df = 1. 



l-be,{L,) 
AC)/.. _ cexiPi - Qfi) 
^1 l-de,(LO 



(4.36) 
(4.37) 



where Li — length(Ji) = P2 — Q^i- Then 



(B,C)= / 4^)(A)AP(A)|5/ii(A)|'dA 



= (ac) 

In particular, if B — C, we get 

(B,B) = |a|^ / - 
Jr 1 

where 6 := |6| e{(p). Also see (4.39). 
Corollary 4.12. Let B = 



(l-6e(ALi)) (l-de(-ALi)) 



2|6|cos(27r(^ + LiA)) + |&|' 



\Shi{\)\'d\. 



\Shi{\)\^d\ 



(4.38) 



(4.39) 



a b 

—b a 



e ;7(2) be as above {b = \b\ e{ip), \a\^ + \b\^ = 1), then 



(PER\Shi{-)\^^ (A) = L\. (See [BJ02].) 



For the Poisson-kernel 



we have: 



n(0 := 



l-\b\' 



l-2|6|cos(27rO + |6r 

1 

{B,B)=Ll Pbi^ + LiX)dX = Li, 



(4.40) 



(4.41) 



(4.42) 



and 



{B,C)^[a-c)L\ 



dX 



(l-6e(ALi))(l-de(-Aii)) 



1 - &d 



^1 



(4.43) 



for all B = 



a b 
—b a 



and C = 



c d 
—d c 



e u{2). 
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Proof. See Example 4.11. 

The justification for the identity (4.40) is from wavelet theory. Indeed, the assertion (4.40) is 
equivalent with the fact that the Shannon wavelet is an ONB- wavelet in L^(E). The summation 
in (4.44) below is justified since \Shi\ is in L^(K), and as a result PER|5ft-i| is in in any 
period-interval. 

Since the first function inside the integrals in (4.46) and (4.39) is periodic with period 1/Li, we 
introduce a periodized version of the function as follows 

2 



PER|5/ii|^(A) 



Shi A 



Lj, (See [BJ02].) 



Note 



and as a result, we get 



PER|S'/ii|^(A+ -5-) =PER|5/ii|^(A), VA e 
Ll 



(PER\Shin (A) 



{B,C) = {ac) 

= (ac) Ll 



dX 



(l-6e(ALi))(l-de(-ALi)) 
fTk 1 

(1 - be{XLi)) {1 -de{-XLi)) 
{ac)Ll / ' ^ 6™d«e((m-n)ALi)dA 



dX 



(4.44) 
(4.45) 



(ac)L? £(M)" dX 



1-bd 



Ll. 



Eq (4.42) follows from this. 

A direct verification of (4.42) is given as follows: 



{B,B) = \a 

= \a\'Ll 
^Ll 



(PERlShin (A) 



l-2|6|cos(27r((p + LiA)) + |&|^ 
1 



.dX 



rdX 



l-2|6|cos(27r((p + iiA)) + |5|' 

2 

2*^^ = Ll. 



(4.46) 



(4.47) 



/o l-2|6|cos(27r((^ + LiA)) + |6|' 

The conclusion (4.42) follows from (4.47), (4.41), and the normalization property of the Poisson 
kernel. ■ 

Remark 4.13. Eqs (4.40) and (4.44) can be verified as follows (also see [BJ02]): Let g := xJn and 

g{x) :— g{~x), so that 



{g*~gr{X)^\ShiWX) 



sin(7rALi) 



ttA 
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It follows that 



Jo ^ ^ 7r 




n^Q. 



Hence, |^PER |S'/ij|^j (A), as an j^-periodic function, has Fourier series 

(PER|5/ii|') (A)=co = L?. 

4.2. Unitary Dilations. For the definitions of key notions, unitary dilation for semigroups of con- 
tractions, minimal unitary dilation, uniqueness up to unitary equivalence, and intertwining operators, 
the reader is referred to e.g., [APOl, J0r81, LP68, Vas07], but the literature in the subject is extensive. 

For each of the selfadjoint extension operators Pg, we now compute three associated scattering 
semigroups. There is one for each of the infinite half-lines contained in 17, and one for the union of 
the bounded components in fi. 

Now, this means that when B is fixed in [/(n) then the unitary one-parameter group Usit) will be 
a unitary dilation of each of the three semigroups, one for the incoming subspace in L^(ri), one for 
the outgoing, and a third for bounded components. This fact yields additional scattering theoretic 
information for our problem. 

Remark AAA (Unitary dilations). Now, as before pick n > 2, and fix the open set f2 with interval 
endpoints as in section 3.1 (see Figures 3.1 and 2.4.) When a boundary condition B is fixed, we get 
an associated unitary one parameter group Usit). We claim that this is, at the same time, a unitary 
dilation of three different semigroups of contractive operators (called "contraction semigroups.") 

For every B G U{n) we have a strongly continuous unitary representation Ub{-) of (M, -t-), i.e., of 
the additive group of M. When t g M is fixed, UB{t) is a unitary operator acting in L?{^). 

Associated with {U B{t) , L'^ (Vl)) , one has three contraction semigroups. Each of the three semi- 
groups is the result of cutting down Usit), i > 0, with three separate projections. 

First, let P_ be the projection onto L^{J-), see Figure 2.4, let P+ be the projection onto L^(J_(-), 
and finally Pm denotes the projection onto "the rest," i.e., onto of the union of the n — 1 finite 
intervals Jj, see Fig 2.4. 

The three semigroups are now as follows: 

(1) Z^{t) := P^UB{t)P-, t > 0; its infinitesimal generator has von Neumann indices (0, 1). 

(2) Z+{t) := P^UB{t)P+, t > 0; its infinitesimal generator has von Neumann indices (1,0) 

(3) Zjn{t) PmUB{t)Pm, t > 0; its infinitesimal generator is maximal dissipative, see [JM80, 
LP68]. 

Moreover one easily checks ([Kos09, JM80, LP68]) that the first is a semigroup of co-isometries, the 
second a semigroup of isometrics, and the third one Z„i{t) is a semigroup of contraction operators; 
often called the Lax-Phillips scattering semigroup. 
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Figure 4.2. The action of the two semigroups, Z-{t) acting on L^(— cjo, /3i), and 
Z+{t) acting on L^(q!„, oo), i > fixed. The semigroup Z+{t) is called a semigroup 
of shifts: Each operator Z+{t) shifts functions to the right. A given function h+ on 
the halfline (q;„,oo) is shifted like a Xerox copy. Z^{t) moves it t units to the right. 
The norm is preserved: It is a semigroup of isometrics. By contrast, moving 
instead to the left in (q;„, oo), we get a semigroup of co-isometries; hence truncation, 
and the norm is decreased. By symmetry, this duality works in the reverse on 
the LHS half- line (— cx),/3i). In general the adjoint of a semigroup of isometrics is a 
semigroup of co-isometries. 



Lemma 4.15. IfY± denote the respective infinitesimal generators of the two contraction semigroups 
{Z±(t)}tem^ in (1) and (2) above, then for their respective dense domains we have: 

&{Y+) ^{f;f and f e L^K, oo), /(a„) = 0} (4.48) 

dense in — L'^(J^), J_|_ — (a„,oo); and 

^{Y-) ^{f;f and f e L2(oo, /3i)} (4.49) 

dense in — L?{^JJ), J_ — (— oo,/3i). 

As before the interval systems a = (a^), and (3 = {j3i) are specified as in Figure 3.4-. 

Proof See e.g., [LP68]. ■ 

Remark 4.16. Using [Kos09], one further checks that ([/^(i), L^(f2)) serves as a unitary dilation of 
all three semigroups. Recalling [Kos09] finally, that, in general, unitary dilations may or may not be 
minimal. As an application of Theorem 4.1, it follows that ([/^(f), L^(f7)) is a minimal dilation of 
Z-{t), t > 0, if and only if B is non-degenerate. In this case, it is also a minimal unitary dilation for 
the semigroup on the right Z^{t). 

But minimal unitary dilations are unique up to unitary equivalence [Kos09]. As a result, we get a 
system of unitary intertwining operators. Below we outline formulas for these intertwining operators, 
and their relevance for scattering theory, and for bound-states. 

Corollary 4.17. If B £ U{n) is non-degenerate, then the scattering coefficients A)^ in (4.2) satisfy 
A^'^'iX) > for all A G M. Moreover if Bi,B2 G U{n) are both non-degenerate, then there is a 
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unitary intertwining operator W in L'^{fl) subject to the following two conditions: 

iwh = h, yh e L2(-oo,/3i) 

\WUBAt)h - UB,{t)h, Vt e e L2(_oo,/3i). 

In Fourier domain, it is determined as follows: 



(4.50) 



( 




9] (^) 



(4.51) 



for 1 < J < — 1, ,9 G L^(M), X G fl; where — ^* g denotes the inverse Fourier transform in 



In this section we make a comparison between families of self adjoint extensions that have purely 
continuous spectrum, and the cases with embedded point-spectrum. We outline detailed scattering 
properties; and in particular, we give examples of non-periodic periodic spectrum; see Theorem 5.6, 
and the caption in Figure 5.3. 

The following features from quantum theory are reflected in property certain of our operators Pg 
which allow decomposition (degeneracy); see e.g.. Theorem 3.20 above, and Figure 2.4; as well as 
details in the section below. States in quantum mechanics are represented by wave functions, and they 
in turn by vectors (of unit-norm) in Hilbert space. In the paragraphs below, we will use wave-particle 
duality (from quantum theory) without further discussion; hence referring on occasion to particles 
as opposed to wave functions. Bound-states are states that satisfy some additional confinement 
property. Consider now a quantum system where particles (waves) are subject to confinement, by a 
potential, or by a spatial region, e.g., a box, or an interval. These particles then have a tendency to 
remain localized in one or more regions of space. In our present analysis we will identify this case 
by a discrete set of spectral-points, embedded in a continuum spectrum (embedded point-spectrum); 
hence the presence of eigenvectors for the operator Pb under consideration. The spectrum, for us, 
will refer to Pb, one in a family of selfadjoint operators (quantum mechanical observables, such as 
momentum, or position). 

In quantum systems (where the number of particles is conserved), a bound-state is a unit-norm 
vector in a Hilbert space which is also an eigenvector. They may result from two or more particles 
whose interaction energy is less than the total energy of each separate particle. Hence these particles 
cannot be separated unless energy is spent . The mathematical consequence is that the corresponding 
energy spectrum for a bound-state is discrete, and in our case, embedded in continuous spectrum. 
Bound-states may be stable or unstable, and this distinction will be illustrated for our model below. 
Positive interaction energy for bound-states corresponds to "energy barriers," and a fraction of the 
states will tunnel through the barriers, and eventually decay. Stable bound-states are associated to, 
among other things, stationary wave-functions, and they may show up as a poles in a scattering 
matrix (details below.) 

Fix n = 3. Recall that £7 — IJj^Q Ji, where Ji — (ai, /32), J2 — ft), J- = Jo = (—00, /3i), and 
J+ — J3 — (as, 00). For B G U{3), the generalized eigenfunction is 



(M) 



5. Degenerate Cases 




(5.1) 
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^0 ^ 




( A, 




\-\ 


A2 


Ai ) 




\ ^3 



with Xi :— xji, « = 0, 1,2,3. The coefficients (Ai)^^^ satisfies the boundary condition 
as in (3.23). 

5.1. Case 1: With [/^(t) indecomposable. Let 

B = 



(5.2) 




et/(3) 



where \a\^ + = 1. Then 



B' 



b 
a 



, w 



and B'*B' 



1 




I -P.. 



Note that = 1. 

Standing assumption < |6| < 1. The notation used in the example is the one introduced above. 

Conclusions: After a computation we arrive at a closed-form formula for all four generalized 
eigenfunction (GEF) coefficients Ai, with the local index i from to 3; and, as a result, a closed-form 
formula for the GEFs ipx in Theorem 2.4; see also (4.10). 

To ensure that the measure in the spectral representation ub is Lebesgue measure, we pick Aq = 1. 
Some noteworthy properties of the coefficients: The coefficient Ai for the first of the finite intervals 
inside f2, carries more information than the remaining three coefficients. Studying transformation of 
states in L^(r2) under unitary one-parameter group Usit), with t increasing, we note that the last 
GEF-coefficient ^3(A) measures transition into the infinite half-line to the right. It turns out that the 
last coefficient, A^ is just a phase factor times the unitary scattering operator A2. All coefficients, 
phase factors, and time-delay depend on the respective lengths of the finite intervals in ^l, as well as 
the lengths of the gaps between them. 

The A2 function is a scattering operator (see [JPTllb]) adjusted both with a phase factor and 
an additive time-delay. Hence, three of the four GEF-coefficient have modulus 1, i.e, \Ai (A)| — 1 , 
for i — 0,2 and 3. The coefficient |^i| carries a probabilistic interpretation. It is a scaled Poisson 
kernel, with the scaling depending on a two-state distribution \a\ + \b\ = 1, where a and h are 
complex, the SU{2) entries from B. 

As a result, in the spectral decomposition (Theorem 4.1), we get local densities — 1, except at 
one place, for first of the finite intervals Ji, where the distribution density is \Ai\ . So by contrast 
to the case n — 2 [JPTllb], in the present model, we do not have Poisson uniformly contributing to 
as- The spectrum of Usit) is pure Lebesgue spectrum, with no embedded point-spectrum. And the 
global Hilbert space L'^{n) does not decompose. 

The boundary conditions are as follows: 





In details, we get the following transformations: 
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/(ai) =a/(/3i) + 6/(^2) 
/(a2) =-6/(/3i)+a/(/32) 




V 



Figure 5.1. Transition between intervals in Q, 



Setting 



so that 



and we have 



Le., 



I aex{Pi-a-C) bex{/32-ai) 
= -bex{l3i-a2) aex{P2-a2) 

V 



exiPs - as) / 





a:i) 


aex{P2 - 


0:2) 


^0 \ 


(A, 


A, = 


A2 


A2 J 


\ As 



aexiPl 
-bexiPi- 



' o-i)Ao 
0:2)^0 - 



-bexW2 
aexih 
ex{Pz 



ai)Ai 
■ a2)Ai 
as)A2 



Ai 
A2 
^3. 



The characteristic polynomial of B'^ ^(A) is 

X — b ex{l32 - ai) 
-aex{l32- Oi2) x 



det 



x{x-bex{h)) 



with roots x = 0, x = bex{li) 7^ 1. But since |6| < 1 we get 1 ^ sp{B'^^^{X)), VA e K, hence Apt 
and {I2 — Ba^p{Xj)~^ is well-defined for all A G M. 
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Note that 







1 -hex{li2 - ai) \ aexih - "2) 1 ~bex{f32 - ai) 

Setting = 1, it follows that 



-bexiPi - a2) J 

aex{l3i - ai) 



A 



3 



l-bex{l32- ai) ^ ex{Pi + ^2 - ai - a2) - bex{Pi - a2) J 

Finally, 

1 // \ / aexWi-ai) \\ 

1 - bexW2 - ai) \ V eA(/33 -as) J'^ g^(^^ + /^^ _ - ^2) - SeA(/i) / / 
= exik - 0:3)^2 
We summarize the results in the lemma below: 
Lemma 5.1. The solution to (5.3) is given by 

Ao = 1 (5.4) 

_ aexWi-ai) , . 

l-bex(A-a,) ^^-^^ 

. _ exiPi - a2) {ex{l32 - ai) -h) 
1 - bex(p2 - ai) 

As^ex{lis~az)A2 (5.7) 
Lemma 5.2. Setting b = \b\ e{Lp), (p eR, then 

2- (5-8) 

1 - 2|6|cos(27r(v3 + /iA)) + |6|^ 

Hence \Ai\'^ is the Poisson kernel with parameter b. 

Proof. This follows from (5.5). I 

5.2. Case 2: With UB{t) decomposable. Let 

a b 

B = \ -b a \ e U{3) (5.9) 

1 
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where 



1: and u = w 



B' 






(5.10) 



and c — 1. 



Summary of conclusions in the Example. Standing assumption < |5| < 1. 

The notation used in the example is as before, but the element B in [/(3) is now different. We 
also fix a system ot and /3 of interval endpoints, subject to the standard position, see (2.1) through 
(2.3) in Section 2.1. 

As always, the conclusions will depend on both B and the prescribed pair a. and (3: Again, we 
arrive at a closed-form formula for the generalized eigenfunction (GEF) see Theorem 2.4 and 
(4.10). But this time, we get embedded point-spectrum in the continuum (boundstates in physics 
lingo.) 

The discrete set Apj making up the point-spectrum depends on both the lengths of the two finite 
intervals Ji and J2, as well as on the gap between them, and the gaps to the infinite half-lines. 

As before, to get the continuous part of cjb to be Lebesgue measure on M, we pick Ao — 1. 

Studying transformation of states in L^(i7) under unitary one-parameter group Usit), with t 
increasing, we note that incoming states from the infinite half-line to the left turn into boundstates. 
But the action of [/^(t) on the global Hilbert space _L^(il) now decomposes as an orthogonal sum of 
continuous states, and boundstates. 

As a result, in the spectral decomposition (Theorem 4.1), we get local densities — 1, for the 
continuous part, and a set of Dirac-combs for the discrete part. But by contrast to the case n = 2 
[JPTllb], in the present model, we get non-periodic Dirac-combs. The spectrum of Usit) is a mix 
of Lebesgue spectrum and embedded point-spectrum. 

The boundary condition takes the form 

/(«i) 

/(«2) I = 

/(as) 




so that 




= a/(/3i) + &/(/32) 
= /(/3i). 



(5.11) 



See the second line in Fig 5.2 (and also Fig. 2.4) for a geometric representation of the last equation 
/(as) = /(/3i) in the system (5.11) of boundary conditions. 
In this case, 



Ba,pW — 



( ° 



V 1 



aex{l32-ai) bex{l33-ai) 
~-be\{l32 - a2) ae\{l3z-a2) 




where 



~b ex{fi2 - 0-2) aex{l33 




e SU{2), VA e 
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L^in) = L2(J_ U J+) ® L^{Ji U J2) 

Figure 5.2. Transition between intervals in fi. 

Now, the boundary condition (5.2) becomes 

aex{f32 - ai)^i ^bexiP^ - ai)A2 = Ai 
-bex{fi2 - a2)Ai + aeA(/?3 - a2)^2 = ^2 

We set ^o(A) = 1 for all A € M. 
As a result, we see that the vector 

A(\\ — I 

A2iX) 

must be an eigenvector of B'^ ^(A) for A to be in the spectrum of Pb, or equivalently for the unitary 
one-parameter group J7s(i). 
As a result we get 

Apt = {\eR; detih - B'^JX))-' = 0}. (5.13) 
Hence, as the interval endpoints a ~ (ai) and /3 ~ {Pi) fixed, the set Apt results as the solution 
manifold for 

/ 1 - aeA(/32 - "i) ~h ex{l3:i ~ ai) \ 
det ( _ =0. (5.14) 

\ bex{(i2-a2) I - aex{Ps. - a2) ) 

Notice (5.14) is independent of /3i and as. 



^(A)=(' tS'! ) (5.12) 



Example 5.3. Let 



a = {1,2, 3 + 93}, > 
3 

^ = {0, 2 '3}, and 

a = b = 

V2 



Substitute into (5.14) 



det I ^ ' ^ =i_ 1 (lA)-i=e(A) + e(^A). 
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As a result 



z(A) 



e(-A) — cos(7rA) + i sin(7rA) 



must satisfy the following cubic equation 



1 - - + = (1 + Z){7? - (1 + 



Hence, 



or 



(1 



V2 



z = -\ 



U + 1 = 0' 



V2 

A e 1 + 2Z; 



z + l) =0. 



(5.15) 



Note \z± \ — 1. Let A± be such that z± — e(|A±), with A± e M. We conclude that 

Apt = (1 + 2Z) U {A± + ^Z}. 

Remark 5.4. In any example with as in Figure 5.2, i.e., when is the complement of three finite 
closed intervals, O will have two bounded components, i.e., open intervals J^, i = 1,2; and two 
unbounded. If further Usit) is assumed decomposable, there will be one summand U'"^''*B{t) acting 
on L^{Ji U J2) of the union of the two intervals J^. 

Example 5.3 produces one particular configuration for this possibility, so a computation of the 
spectrum of U^'''^* B{t) when there are boundstates. In an earlier paper [JPTlla] we found all the 
configurations for spectrum for each one of the possible momentum operators in LP' of the union any 
pair of finite open intervals Ji. 

This in turn is a question of interest both for the study of both quantum systems, and of spectral 
pairs, see e.g., [Fug74, JP98, JP99, DJll, LabOl]. 

Remark 5.5. More generally, if 



( 





g 









V ••• 


1 



e U{n) 



with g e SU (n — 1) , then the corresponding unitary one-parameter group Ubi (t) does not decompose. 
See case 1 of section 5.1. 



On the other hand, for 



( 







V 1 



g 







e V{n) 



/ 



where g G SU{n — 1) as before, the unitary group UB2{'t) decomposes. See case 2 of section 5.2. 
Set 

/ ••• 1 \ 



1 



vo 







1 oJ 
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then (see (3.47)) 



/ 





g 



5* = 



V 1 



y 



g 



V 







1 / 



Theorem 5.6. Let Li = length{ Ji), i = 1,2, be the lengths of the two bounded intervals Ji and J2. 
Then, 



D{X,a,Li,L2) :=dct(B^^^(A)) 

= 1 + ex{Li + L2) - aex{Li) - aex{L2); 



(5.16) 



where det(i3^ pi^)) defined in (5. 14-)- The solution manifold Apt (5.13), i.e., the embedded point- 
spectrum, is the set of zeros 

Z{a,Li,L2) :={AeIR; D{\,a,Li,L2) =0} (5.17) 
of the exponential polynomial in (5.16). Moreover, setting a :~ we{ipQ), < w < I, eq (5.16) is 



equivalent to 



e{XL2 - <y3o) = 



1 — u; e(ALi + fo) 
w — e{XLi + ipo) 



(5.18) 



Proof. Eq. (5.16) follows from a direct computation. As noted in [JPTlla], both sides of (5.18) can 
be interpreted as periodic motions on the torus T^: 

(i) The LHS is a uniform motion with constant velocity; 

(ii) The Mobius transformation on the RHS has the form e'^^^^, where 



rdu. 



(5.19) 



2 ' 27r Jo 1 2 Jq 1 ~ 2w cos(2ttu) + 

The solution to (5.18) is obtained at the intersection of the two motions. In particular, the solution 
(point-spectrum) is periodic if and only if L2/L1 is rational. See Figure 5.3. I 

Remark 5.7. For all A G Apt = Z{a, Li, L2), the vector of coefficients (5.12) satisfies 

Ai (A) = '"-J-— e (A (/?3 - ai)) A2 (A) 

1 — a e (XLi) 



A2 (A) 



e(A(/?2-a2))Ai (A). 



1 — a e (AL2) 

Note that \Ai (A)| = |^2 (A)|. As shown in [JPTlla], the pairs of sets 

(JiU J2,Apt) 

forms a spectral pair if and only if 

A, (A) = A2 (A), VAeApt; 

and when (5.21) holds, Ji U J2 is said to be a spectral set. 

Indeed, the union of [1, |] and [2, 3] in Example 5.3 is not a spectral set; one way to see that is to 
notice that it is not a tile for the real line under translations: you cannot fill the gap [|, 2]. For more 
details, see [JPTlla]. 



(5.20) 
(5.21) 
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Figure 5.3. Dense orbits: The spectrum as a random-number generator [JPTlla]. 
This is obtained in a number of steps: Step 1. Consider the curves for argument 
function from the right-hand side of (5.18) and the hues with slope L2 from the left- 
hand side of (5.18). Assuming L2/L1 is irrational. Step 2. Identify the asymptotes. 
Let Si, i G Z, he the i*'' interval between neighboring branch cuts; fix Sq to be 
the closest to 0. Note that all the intervals have the same fixed unit-length, and 
they extend both to the left and to the right of Sq. Step 3. The embedded point- 
spectrum Apt is discrete and infinite, it intersects all these intervals between branch 
cuts as Ai := Apt n Si. Now, translate all of these finite intersections down to 5*0. 
Conclusion: Since the line-slope L2/L1 is irrational, the set Uigz(Ai — i) is dense in 
So. 



5.3. Other Examples. 



Example 5.8. For n = 3, let 
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and so B'^ 
equation 

we get 



eA(ii) 







where Li ~ length(Ji), i = 1,2, as before. Solving the 



-exiL^) 
dct (/2 - B'^j,{X)) = (1 - GA (ii)) (1 + ex {L2)) - 



Apt = 



1 



Z U 



-Z 



As shown in Figure 5.4, the Lebesgue spectrum arises from lumping together i^(J_) and L^(J+); 
and the embedded point-spectrum Apt accounts for the boundstates in i^(Ji) L^{J2). 



Ao 
J-=h 



Ai 
^1 



A2 



»2 



A3 



ySi 0:3 

Figure 5.4. Embedded point spectrum. 
Example 5.9. Let n = 3, and 

1 
S = I 621 623 I e U{3) 
631 633 



assuming I623I ^ 1- Here, B'^^^iX) = [ f'^ b^s^xiL^) 



The determinant criterion 



det 



ex (Li) - 1 







yields 



b23 ex (L2) - 1 

1 „ 



Apt = 



As illustrated in Figure 5.5, we have 

acting on ^^(Ji) ® L2^J_ U J2 U J+). For a detailed analysis of f/|3°"*(t), see [JPTllb]. 




Figure 5.5. Embedded point-spectrum. 
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6. Decomposability 



As we outlined in sections 3 - 5, as B g U{n) varies, the unitary one-parameter groups Usit) act 
in L^(f2). Now the given open subset is a disjoint union of its connected components, i.e., of a 
specific set of intervals. As a result, L^(0) splits up as an orthogonal direct sum of a corresponding 
number of closed subspaces; one L^-space for each of the component intervals. But it is also true 
that that the typical scattering theory for Usit) corresponds to an action in L^(ri) that mixes these 
closed subspaces in Indeed, when B G U{n) is fixed, our results Corollaries 4.6, 4.8, 4.10; 

Remark 4.14, and Figure 4.2 yield formulas for transition probabilities, referring to transition between 
the interval-subspaces, and governing the global behavior of UB(t) as it acts in L?(yi). The term 
"decomposability" in the title above refers to invariance under UB{t), for all i € M, of some of the 
interval-subspaces in i^(ri); clusters of subspaces. 

In this section it is convenient to use a slightly different labeling of the selfadjoint operators Pg. 



Let il := Ufc^o-^fc where Jq :=]-oo,/3„], Jk ■= [afc,/3fc], fc = 1,2, 



1, and Jn := c>o[. So fl is 



the complement of n intervals: 17 = M \ Ufe=i]/^fe' o^fef- The selfadjoint restriction of P are indexed by 
the unitaries B from ^^(afc) — ?► P{j3k). Identifying the spaces P{ak) and £'^{f3k) with C" we realize 
i? as an n X n matrix. 



J- = Jo 



Pn «i /3i 



0-3 



Pn- 



Figure 6.1. The complement of n bounded intervals in M (n > 2). 



As usual the domain of the maximal operator is the absolutely continuous functions on Q and the 
selfadjoint restrictions Pb are in one-to-one correspondence with the unitaries B. The domain of the 
selfadjoint restriction Pb determined by B is the set of absolutely continuous functions / : — >■ C 
satisfying the set of boundary conditions 



B 









fia2) 













(6.1) 



and Pb/ — j^f- Suppose B is block diagonal, this is 



B = 



Bi 
B2 



= Pi ® P2 



where Pi is a fc x fc matrix and P2 is a (n — fc) x (n — k) matrix. Then Pj, j = 1,2 are unitaries and 
we can write 17 = ili U where 



17i := Ji U J2 U • • • U Jk and 

^^2 '■— Jq U Jk+i U Jk+2 U • • • U J„. 
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Consequently, (il) 
determined by 



LP' (rii)©L^ {^2) and Pb = Pbi ®Pb2 where Pb^ is the momentum operator 



Bi 







r fWi) 1 


/(«2) 




/(/32) 









(6.2) 



and Pb2 is the momentum operator determined by 



Bo 



Hence, if B is block diagonal it is sufficient to study Pb^ and 



/("fc+l) 




" /(/3fc+i) " 


/(afe+2) 




fWk+2) 


/(an) 







(6.3) 



Remark 6.1. A reason for grouping the unbounded intervals this way is that the deficiency indices 
work this way. The restriction of P to each {Jk) , k = l,...,n — 1 and the restriction of P to 
(Jo U Jn) all have deficiency indices (1, 1). Consequently, the restriction of Pb^ to (fli) has 
deficiency indices {k,k) and the restriction of Pb^ to (^2) has deficiency indices {n — k,n — k). 
Furthermore, k = 2 the Pb^ problem is investigated in [JPTlla] and if n~ k — 2 the problem 
is investigated in [JPTllb]. 



Recall, a permutation matrix is an nxn matrix obtained from the identity matrix 
by permuting of the columns of /„. 



diag(l, !,...,!) 



Definition 6.2. We say two unitary matrices A and B are permutation equivalent, if there is a 
permutation matrix S such that B = S*AS. We say a unitary matrix B is decomposable, if B is 
permutation equivalent to a block diagonal matrix and we say B is indecomposable, if B is not 
decomposable. 



Example 6.3. Let n — 4, and 



11 





bi3 








^22 





^24 


31 





&33 








h2 





644 



G C/(4). 



The boundary condition reads 



'11 





^13 







" /(ai) " 




' fWi 





b22 





624 




/(a2) 




fW2 


'31 





^33 







/(as) 




/(A 





^42 





644 _ 




. 1(^4) _ 







Note that i? is permutation equivalent to 



A 



and it follows that the system decouples as shown in Fig. 6.2. 



'11 


bi3 








'31 


bss 














^22 


624 








^42 


644 
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^44 




^42 fon 



Figure 6.2. B is permutation equivalent to A. The system decouples into a direct 
sum of two subsystems: (i) The dashed diagram contains two bounded intervals Ji 
and J3, and this corresponds to [JPTlla]; (ii) The solid diagram consists of one 
bounded component J2 and two unbounded components J±; and it is investigated 
in [JPTllb]. 



Supposing B is permutation equivalent to A we can write (6.1) as 



AS 







\ m) 1 


/(«2) 




m) 


= s 









using that 5" is a permutation this can be written as 



So Suppose B is permutation equivalent to a block diagonal matrix A = © A2, then B is permu- 
tation equivalent to a block diagonal matrix A such that i„ = n. Putting it together we have 

Theorem 6.4. If B is decomposable, then we can write S*BS = BiQ) B2(B ■ ■ -©-Bfc where each Bj is 
indecomposable and S is a permutation. The Pbj , j = 1, 2, . . . , fc — 1 problems only contain bounded 
intervals and the Pb^ problem contains the unbounded intervals and, perhaps some of the bounded 
intervals. 

Since the unbounded intervals are "special" it is useful to write as 

w c 

Lemma 6.5. If B is decomposable, then B' is degenerate, i.e., has an eigenvalue with absolute value 
one. 
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7. ElGENFUNCTIONS 

Fix some unitary matrix B. The generalized eigenfunctions 

ipx{x) := Ao{X)x]^ac.,p,,[{x) + ^ Aj {X)xiaj, 13 jlix) + An{X)xia„,oc.[{x) ex{x) (7.1) 

that satisfy (6.1) for the generalized eigenspace corresponding to A. The coefficient Aj = Aj{X) is 
obtained by solving the differential equation -^ip — 2niip on the interval Jj. Plugging (7.3) into (6.1) 
we see the generalized eigenfunctions are determined by the solutions Aq,Ai,--- ,An to the system 
of n linear equations in n + 1 unknowns: 



B 



Aie{Xai) 
A2e{Xa2) 

Ai-ie(Aa„_i) 
Aae{Xan) 



Aie(A/3i) 
A2e{X/32) 

^„_ie(A/3„_i) 
Ane{X(3n) 



(7.2) 



Let Da ■- diag(e(Aai),e(AQ;2), ■ • • ,e(Aa„)), := diag(e(A/3i), e(A/?2), ■ . ■ , e(A/3„)), and Ba./3 := 
D*^BDa- Then our eigenvector equation can be written as 



Writing C" 



a. 13 



we have the decomposition 
B 



Ai 




Ai 


A2 




A2 


An-l 




An^l 






A„ 



(7.3) 



■>a,l3 



B' u 
w c 



(7.4) 



notation we can write (7.3) as 



u, w 


are 


m 




and B' is 


" B' 


u 




■ A' ' 




■ A' ' 


w 


c 




. ^0 . 




An 



and B' is a (n — 1) x (n — 1) matrix. With this 
where A' ^ [Ai, A2, . . . , A^-i] . 



Theorem 7.1. If \i is in the range of T — B' and rjQ is such that u = (/' — B') rjQ, then the solutions 
to (7. .3) are determined by: 

A' = Aaija + ( 
An = Aq{c + w?7o) + wC, 

where C £ ker (/' — B') and Aq G C. If u is not in range of I' — B' , then the solutions to (7.3) are 
determined by: 

A'^C 



where ( e ker (/' - B') and Aq = 0. 
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Proof. We can write (7.3) as 

B'A' + Aqu = A' 
•wA' + cAo = An- 

The result is immediate from this. ■ 

Neither the theorem not the first corollary require B to be unitary, but the second corollary needs 
u = implies w = 0, which is a consequence of the assumption that B is unitary. 

Corollary 7.2. //I is not an eigenvalue jor B' ^ then the solutions to (7.3) are determined by: 

A' = Ao [I' -B'Y^u 

A^ = A^ (^c + w{I' -B'y'uj . 
In particular, the set of solutions to (7.3) is one dimensional. 

Proof. If 1 is not an eigenvalue for B' , then the kernel of /' — B' equals {0} and the range of /' — B' 
is C"^^, in particular, u is in the range oi I' — B' . ■ 

Corollary 7.3. If u — 0, then the solutions to (7.3) are determined by: 

A'^C 

An = Aqc + wC 

where ( e ker (/' - B') and Aq G C. 

Proof. If u = 0, then u is in the range of /' — B' . Since B is unitary |c| = 1, hence w = 0. ■ 

An immediate consequence of Theorem 3.9 and Corollary 3.28 we have 
Theorem 7.4. // B' is not degenerate, then the spectrum of Pb has uniform multiplicity one. 

Corollary 7.5. Suppose B is decomposable with decomposition ® Bj®Bk in the sense of Theorem 

6.4 and each Bj is not degenerate, then the spectrum of Pb — Pbj ® Pb^ where the spectrum 

of Pbj is a set Aj of simple eigenvalues and Pb^. has spectrum equal to the real line and the spectral 
measure is absolutely continuous with respect to Lebesgue measure. 

In particular, the set of eigenvalues of Pb is Uj=i ^3 ^-^i^ the multiplicity of an eigenvalue A is the 
number of elements in {j = 1, 2, . . . , fc — 1 | A G Aj} . 

8. Scratching the Surface of Infinity 

In this section we consider some cases when the give open set VI has an infinite number of connected 
components. As in the discussion above, we still assume that two of the components are the infinite 
half- lines. Our motivation for studying the infinite case is four- fold: 

One is the study of geometric analysis of Cantor sets; so the infinite case includes a host of 
examples when 17 is the complement in M of one of the Cantor sets studied in earlier recent papers 
[DJ07, DJll, JP98, PWOl]. The other is our interest in boundary value problems when the boundary 
is different from the more traditional choices. And finally, the case when the von Neumann-deficiency 
indices are (00,00) offers new challenges; see e.g., [DS88]; involving now reproducing kernels, and 
more refined spectral theory. 
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Finally we point out how the spectral theoretic conclusions for the infinite case differ from those 
that hold in the finite case (see details above for the finite case.) For example, for finitely many 
intervals (Theorem 3.20) we computed that the Beurling density of embedded point spectrum equals 
the total length of the finite intervals. By contrast, we show below that when D, has an infinite 
number of connected components, there is the possibility of dense point spectrum; see Example 8.5. 

Let Ik = (rkjSk) be a sequence of pairwise disjoint open subintervals of the open interval (0, 1). 
Let 

■DO 

(-cx),0)U(l,c^)U y Ik. 

k=0 

The functions satisfying the eigenfunction equation ^^"^a = X4'\ £^re the functions 

Tpxix) = I A_oo(A)x(-oo,o)(a;) +^oo(A)x(i,oo)(a;) + ^ Ak{X)xi^ix) j exix), 

\ k=0 ) 

where j4_oo, ^oo, and Ak are constants depending on A. Let rg = 1 and sq = 0. 

Example 8.1. An example of this is the complement of the middle thirds Cantor set C. We can 
write the complement of the Cantor set C as 

oo 2^ 

(-oo,0)U(l,oo)U U U (a,^fe,a,-fe+3-(J+i)) 

j=0 fc=l 

ao,i = .1, 
ai,\ = .01, ai,2 = .21 

02,1 = .001, 02 1 = .021, 02,3 = .201, 02,4 = .221 



where in base 3 



and so on. So Oj.fc, fc = 1, . . . , 2^ are the numbers with finite base three expansions of the form 

0.a;iX2 • • • Xj\, xi G {0, 2}. 
In this case the generalized eigenfunctions are 

(oo 2^ \ 
^-oo(A)x(_oo,o)(a;) + ^oo(A)x(i,oo)(a;) + ^^^j,fe(A)x(„^.^„^^+3-o+i))(x) ex{x). 

Consider a selfadjoint restriction Pb of the maximal momentum operator on such that Ak £ i'^ 
and 



BDriX) 



where 

DriX) = diag (e(Aro), e(Ari), e(Ar2), • ■ ■ ) = diag (e(A), e(Ari), e(Ar2), • • • ) 
Ds{X) = diag (e(Aso), e(Asi), e(As2), • • • ) = diag (1, e(Asi), e(As2), ■ ■ • ) 
and B is some unitary on i"^. 



^oo 




^-oo 


Ai 




Ai 


A2 




A2 


As 


As 
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Theorem 8.2. If B = diag (1, 1, . . .) , then the spectrum of Pb is the real line and the embedded 
point spectrum is Ap = UfcLi I^^' '^here £k = Sk — rk is the length of Ik- The multiplicity o/ A G Ap 
equals the cardinality of the set {k \ Xik € Z}. 



Example 8.3. Some examples illustrating this are: 

(1) If ik = 2~^, then Ap = 2Z. Let Zodd be the odd integers. The eigenvalues in 2'^'Zodd have 
multiplicity k and has infinite multiplicity. 

(2) For the complement of the middle thirds Cantor set Ap = 3Z. The eigenvalues that are 
multiples of 3*"' but not of S'^"'"^ have multiplicity 2'^' — 1 and has infinite multiplicity. 

(3) If ik/'l^j is irrational for all j ^ k, then has infinite multiplicity and all other eigenvalues 
have multiplicity one. 

Corollary 8.4. If B = diag (e(0o)i . . .) , then the spectrum of Pb is the real line and the 

embedded point spectrum is Ap = Ufc°=i (f^ "'^ Z^^) ' ''^^''^^^ = Sfe — fk is the length of Ik- The 
multiplicity of X Cz Ap equals the cardinality of the set {k \ X£k ~ Ok G Z}. 

When we have a finite number of intervals the point spectrum has uniform density equal to the sum 
of the lengths of the intervals, see Theorem 3.20. The following example shows that this need not be 
the case for infinitely many intervals. 

Example 8.5. Suppose B = diag (e(6lo), e(6li), . . .) and 4 = 2'^. Then 2^ {Ok + m) = 2^ {0j + n) if 
and only if 2^+'^ {9k — Oj) — 2^^^ {n ~ m) . Hence, if 9k — 9j is not an integer when k ^ j, then each 
eigenvalue has multiplicity one. Note 2^9k is an eigenvalue for each k. Hence, if 2'^9k Xq then Aq 
is a limit point of Ap. Similarly, by a suitable choice of the sequence 9k, we can arrange that Pb has 
dense point spectrum. 

Theorem 8.6. // we write = C ® then B takes the form 



If the spectrum of B' does not intersect the unit circle, then the spectrum Pb is the real line and each 
point in the spectrum has multiplicity one, in particular, the point spectrum is empty. 
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